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ABSTRACT 


The propagation of em waves in solid state plasmas 
made anisotropic by a dc magnetic field becomes of 
interest in the construction of non-reciprocal micro- 
wave devices. This report presents theoretical and 
experimental investigations of this behavior for the 
development of isolators at millimeter and submillimeter 
wavelengths . 

Part I of this report presents a study of a 94 GHz 
v/aveguide isolator, which consists of an annular InSb 
column in a circular waveguide with a longitudinal dc 
magnetic field and is operated at 75°K. The basic 
operation of the isolator is due to exclusion or absorp- 
tion of em vjaves in the semiconductor. In a mode 
matching analysis the characteristic equation is solved 
in order to determine the propagation constants in the 
waveguide with the annular InSb column. The attenuation 
constant was then calculated by matching the fields at 
the ends of the plasma colTim.n. The four gyromagnetic 
modes excited in the plasma column as well as the four 
evanescent modes excited in the empty waveguide are 
taken- into account in the analysis. 

Part II of this report presents a study of reflec- 
tion beam isolators at millim.eter and submillimeter 
wavelengths using the non-reciprocal behavior in solid 



state magnetoplasmas at room temperature . Non- 
reciprocal reflection is observed when the incident 
electric field is polarized in the plane of incidence 
and propagates transversely to a dc magnetic field. 

Two geometries are considered: GaAs-dielectric inter- 

face at 94 GHz; and a free space-dielectric coated 
InSb interface for submillimeter wavelengths (337 ym) 

A physical explanation for the non-reciprocal reflection 
based on the electron orbits in the plasma is presented. 
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Chapter 1 


INTRODUCTION 


The propagation of em waves in solid state plasma 
made anisotropic by a magnetic field becomes of inter- 
est in the construction of non-reciprocal microwave 
devices . For example , this non-reciprocity may be 
expected in a system in which circularly polarized em 
fields propagate along the dc magnetic field. When 
the direction of propagation is reversed, it is necessary 
to reverse the dc magnetic field in order to achieve a 
situation geometrically equivalent to the original. 

Since the procession of the electrons about the dc 
magnetic field has a definite sense, the electric 
current associated with this procession is also 
reversed with the dc magnetic field. Thus, the 
properties of the plasma are altered, and different 
propagation characteristics for the opposing directions 
may be expected. 

Of course, an anisotropy can be introduced by the 

complexity due to the symmetry and band structure of 

crystals, but we will assume that in the absence of a 

magnetic field, the unbounded solid state plasma is 
•> 

isotropic for electromagnetic wave propagation. T#ien 
the plasma is placed in a dc magnetic field, the field 
causes the electrons to orbit perpendicular to the 



2 


field, and the solid state plasma is no longer iso- 
tropic. In this case the propagation characteristics 
of the medium are no longer described by the plasma 
and collision frequencies , but the cyclotron frequency 
also is a controlling parameter. One of the great 
advantages of utilization of semiconductors in the 
development of non-reciprocal devices lies in the 
flexibility of these parameters. 

Many investigations dealing with anisotropic 
materials, particularly ferrites, have been carried 
■out since D. Polder [1949] and C.L. Hogan [1953] 
showed possibilities of practical non- reciprocal 
microwave devices using the tensor properties of the 
permeability. Recently the solid state magnetoplasma 
has been used as a basis of non- reciprocal microwave 
devices using the tensor properties of the dielectric 
constant.. These non- reciprocal devices can be 
classified into Faraday rotation, field displacement, 
plasma resonance and reflection beam devices. In the 
next few paragraphs various types of microwave 
isolators using ferrites or semiconductors are reviev/ed 

A microwave ferrite isolator that uses Fa'raday 
rotation was developed by B. Lax and K.J, Button [1955] 
Similar devices using solid state plasmas have also 
been' built [Libchaber and Veilex, '19 62; Kuno and 
Hershberger, 1967] . The basis of operation of these 
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devices is that right and left circularly polarized 
(henceforth abbreviated cp) waves travel with different 
phase velocities through the plasma along the longi- 
tudinal dc magnetic field. A plane wave, which can be 
regarded as the sum of right and left cp waves, will 
therefore undergo Faraday rotation. 

The field displacement isolator using a ferrite 
was first discussed by A.G. Fox, S.E. Miller and 
M.T. Weiss [1955] . A solid state plasma in the 
presence of a transverse dc magnetic field was used for 
a field displacement waveguide isolator by M, Toda 
[1964] , and R. Hirota and K. Suzuki [1966] - When a 
transverse magnetic field is applied to a plasma in a 
waveguide, a Hall current is produced by the dc 
magnetic field and microwave electric fields . The 
surface charge on the plasma and the rf magnetic field 
resulting from the Hall current couple with em waves . 
This coupling makes em waves excited strongly on one 
side of the waveguide and weakly on the other side. 
Therefore reversing the direction of propagation inter- 
changes the waveguide side on which the microwave 
power flow occurs. A resistance strip attached to the 
side of the plasma attenuates the reverse wave, and 
hardly affects the forward wave. 

B.J. Duncan and L. Swern [1957] developed a plasma 
resonance isolator using a ferrite in a coaxial 
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waveguide. A similar isolator in which a solid state 
plasma partially fills a circular waveguide was develop- 
ed, by B.R. McLeod and W.G. May [1968,1971] in this 
laboratory. The plasma can be made to exhibit a 
nearly zero permittivity to the left cp wave and non- 

ii 

zero permittivity to the right cp wave. The operation 
of the isolator is due to the exclusion or absorption 
of the microwave signal , depending on the direction of 
rotation of cp waves with respect to the longitudinal 
magnetic field. It was also found that mode coupling 
improved the performance of this isolator [McLeod, 

1968; McLeod and May, 1971] . 

Reflection devices in which the desired em waves 
do not, propagate through or along, the solid state 
magnetoplasma but are reflected off its surface have 
been developed in this laboratory [Seaman, 1969; Kanda, 
et al., 1970]. Non-reciprocal phenomena are observed 
when, the incident electric field vector is polarized in 
the plane of incidence and its propagation is perpendi- 
cular to a dc magnetic field which is parallel to the 
surface of the plasma. The useful frequency range of 
this reflection device is in the far infrared. 

Theoretical and experimental investigations of non- 
reciprocal effects in solid state magnetoplasmas for the 
development of isolators for ,use at millimeter and sub- 
millimeter wavelengths are presented in this thesis. 



The purpose of Part I is to present a study of semi- 
conductor waveguide isolators at 94 GHz. These milli- 
meter waveguide isolators consist of circular wave- 
guide with an annular coliimn of a semiconductor plasma 
in a longitudinal dc magnetic field. The basic opera- 
tion of the isolator is due to exclusion or absorption 
of em waves in the semiconductor. The purpose of Part XI 
is to discuss some anisotropic effects which may be 
used for reflection beam isolators at millimeter and 
submillimeter wavelengths . Non— reciprocal reflection 
from solid state magnetoplasmas can be observed when 
the incident electric field is polarized in the plane 
of incidence and propagates transversely to a dc 
magnetic field. 



Chapter 2 


PROPAGATION OF EM WAVES IN A MAGNETOPLASMA 

In terms of the electronic energy bands of semi- 
conductors, the effect of the carriers may be intraband 
or interband in nature. The intraband effects will be 
simply referred to as free carrier effects. On the 
other hand, the interband effects involve two energy 
bands . Since the photon energies at millimeter and 
submillimeter wavelengths are significantly smaller 
than the energy of the nearest neighboring band 
(typically 0.1~2 eV or more), the interband effects in 
semiconductors can usually be neglected. Moreover with 
sufficiently small dc magnetic fields (-.5 kG) used for 
94 GHz waveguide isolator experiments at 75 °K, the 
difference of the energies between successive Landau 
levels in semiconductors, ftco , is about 46°K in 
temperature units and also that between Zeeman splitt- 
ings, gBB, is about 23®K in temperature units. 

Similarly IR reflection beam isolator experiments at 
room temperature (-300 °K) require the dc magnetic field 
of 15 kG, and is about 138®K and gSB is about 69 °K 
in temperature units . Therefore in our experiments 
throughout the thesis, where kT>>'hu^ and gBB, the 
quantimn effects due to Landau levels and Zeeman split- 
tings can be neglected. 
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The non- reciprocal devices presented in this 
thesis depend in principle on the numerical details 
of the boundary conditions or more generally, of the 
geometrical configuration. In fact guided microwave 
propagation through a plasma waveguide is quite 
different from the plane wave propagation in unbounded 
plasmas . Here wave propagation depends strongly upon 
the waveguide geometries, i.e., the guide wavelength, 
the waveguide cutoff condition and the waveguide modes. 
Therefore the actual performance of the non-reciprocal 
devices could be predicted only as the result of 
lengthy computation. However, since the basic 
properties of solid state plasmas with a dc magnetic 
field in the vzavcguide remain largely unaffected, it 
appears advisable to review in broad terms the 
characteristics of the wave propagation which lead to 
the discovery of the conditions of non-reciprocity for 
em wave propagation in a magnetoplasma. 

2-1 Wave Propagation in Infinite Solid State Magneto- 
plasma 

A convenient and useful stratagem for dealing with 
electromagnetic waves in a solid state plasma in a 
magnetic field is to characterize the medium as a 
complex tensor dielectric constant. This is calculated 
for the simplest case where a mass is isotropic, and a 
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collision time is isotropic and independent of energy. 
The case of energy dependent’ collision time is discussed 
in Appendix A. Anisotropy of the conductivity is 
expected from the fact that the charged particles' are 
forced to rotate about the magnetic field lines with 
the cyclotron frequency, and the electric field of 

the wave develops a component perpendicular to the 
motion. 

In the case of a single type of charge carrier 
and intraband effects only, the equation of motion is 
written in the form: 

m* || = , (2-1) 


including the Lorentz and Langevin forces . Maxwell ' s 
equations are given by 


Vx3 = 


ll- 


Vx* . - || 


( 2 - 2 ) 

( 2 - 3 ) 


It is assumed that the wave varies sinusoidally in time 

- 5 - 

and has propagation constant P, i.e., exp ( jut-P *r) . 

Then making use of the relation, J = nev = aE, 

VxH = jwe^E + OE = jws (K^T + -.-2 )E = jtoe 1 e . (2-4) 

L o L nme o 

o 

We linearize by including in B only the applied dc 
magnetic field and neglecting the rf magnetic field. 



without loss of generality we adopt the Cartesian 
coordinate system where the applied dc magnetic field 
is along the positive z-axis. The tensor dielectric 
constant K then becomes 


-K, 0 


K = 


K. 



0 


0 


0 



(2-5) 


where 


K, = 




' L = K - n (jmT+l)T , 

L jms L m ^ \ 2 , , ^2^ 

O (j(OT+l) +(o)^t) 


a. 


K„ = 


. 2 : 
TO) U) T 

P r C 


gwe 


o 


0) 


[• 


( jWT+l) ( w^t) ^ 


(2-6) 
] , (2-7) 


^11 T 

K, I = V + Li = P_ [ . ] 

j I , L “ gwT+l 


( 2 - 8 ) 


Definitions of all symbols, x and etc. are 

given in the List of Symbols. The relation between the 
propagation constant F and the tensor dielectric 
constant K is obtained from Eq. (2-3) and (2-4) . The 
result is given in the standard ■ textbook (Allis, et al. 
1963) : 

rx(fxE) - k^ic-E = 0 

o 


(2-9) 
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In the Cartesian coordinate system where the dc 

magnetic field lies in the z-direction, and the 

propagation constant F is in the x-z plane with the 

angle 6 from z-axis as shown in Fig. 2-1, one can 

obtain three homogeneous linear equations in E , E and 

^ x' y 


2 2 2 
k^Kj+r cos e 

o 1 


-k^K 
o X 


2 . 

■r sin0cos6 




2 

k K 
o X 


2 2 
k K,+r 

o 1 




-r sin0cos0 


k^K| |+F^sin0 


E 


( 2 - 10 ) 


A non-zero solution exists only when the determinant of 
the coefficients vanishes. Then the dispersion relation 
becomes a biquadratic equation, which is known as the 
Appleton-Hartree formula: 


where 


and 


4 

AF 

+ BF^ + C 

= 0 , 

(2-11) 

A = 

Kj^sin^ 0+K j 

2 

1 cos 0 , 

(2-12) 

B = 

k^K.K sin^ 
o + - 

0+k^K , [K I (l+cos^0) , 

° II 1 

(2-13) 

c - 

k^K.K K, 1 


(2-14) 






= Kj^±jK, . 


(2-15) 
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Since most experiments are performed with, the propaga- 
tion vector either parallel or perpendicular to the dc 
magnetic fields it is most instructive to consider these 
two cases individually and to solve the dispersion 
equation . 

2~2 Longitudinal Propagation (6=0) 

For the longitudinal case, the propagation constant 
r is parallel to the dc magnetic field along the z 
direction. By setting 0=0 in Eq. (2-10) the secular 
equation becomes 


,2., 2 ,2 



k K,+r . -k K 0 


E 

o ox 


X 

, 2 ■ ,2 2 



k^K^ k K|+r 0 


E 

ox ■ o i 


Y 

0 0 k^K 1 1 


E 

o 1 1 


z 


(2-16). 

Then one finds that E^=0, i.e., the propagating X'^rave is 
TEM. By defining the right (+) and left (-) cp waves, 

E^ = E^ljEy (2-17) 

it is possible to decouple the above equations , and the 
propagation constants for the right and left cp waves 
are given by 
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IK. 


• r 

w 1+jx (a}±(0 ) 


}] 


(2-18) 


■ A relatively simple situation occurs when the 
plasma is lossless^ i.e., and the wave frequency 

0 ) is much smaller than the cyclotron frequency 
Then Eq. (2-18) becomes 


2 2 

rf = -k IK- ± + -1 

± o L uoj 

c w 


(2-19) 


In this case the cyclotron orbits are almost closed 


(w T>>1 ) f and the carriers are essentially tied to 

O X 


magnetic lines of force and must move with them. The 
amplitude of oscillation induced by the electric field 
E of the wave in a direction at right angles to and 
E, the Hall direction, is the same for all carriers 
and is independent of mass in a first order approKima-- 
tion. W^en the plasma is not compensated, the nature 
of the wave is dominated by the Hall currents . This is 
called a Helicon wave. When the plasma is compensated, 
the electric current in the Hall direction (the Hall 
current) is very small because equal densities of 
positive and negative charge oscillate in phase with 
equal amplitudes. Assuming that electron-hole inter- 
action is unimportant, one can add their contributions 
independently term by term in the dielectric tensor in 
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the calculation of the propagation constant. This is 
called an Alfv4n wave. 

The semiconductor materials used typically as the 
solid state plasmas are InSb, GaAs , Si and Ge. The 
electron mobility of these materials is typically 
between one and two orders of magnitude greater than 
that of holes. Thus usually only contributions due 
to electrons need be considered in order to investi- 
gate the propagation characteristics in solid state 
plasmas . 

2-3 Transverse Propagation (0=90°) 

In the case of the transverse propagation, the 
propagation constant T is along the x direction, 
perpendicular to the dc magnetic field along the 
z-axis as before Fig. [2-1] . Then three homogeneous 
linear equations in E^, E^, and E^ are obtained by 
setting 6=90° in Eq. (2-10) . 


k K [ 

-k^K 

0 


E 


O X 



X 

k^K 

2 2 
k K|+r 

0 


E 

O X 




y 

0 

0 

2 2 
k K, i+r 


E 



o 1 1 


z 


It is noted that, for the case of the transverse 



15 


propagation f there is a component of the electric 
field along the direction of the propagation by 
virtue of the Hall current. Because of this longitu- 
dinal component of the electric field, the existence 
of space charge may be expected. The dielectric 

relaxation time £/cr, with which any space charge 

-13 -14 

disappears, is order of 10 to 10 sec for the 

semiconductor materials used in the experiments 

described herein. Since we are dealing with phenomena 
. -11 -12 

of period 10 to 10 sec, we can neglect the space 
charge. This is also confirmed from the fact that the 
space charge calculated from the longitudinal component 
of electric field through the Poisson's equation is 
indeed negligible, and hence the charge neutrality 
condition is valid. 

The propagation constant for the case of the 
transverse propagation with the incident em electric 
field parallel to the dc magnetic field is given by 


r 


2 





2 

w 

^ — L_] 

W jWT+1 


.( 2 - 21 ) 


This mode propagates independently of the dc magnetic 
field. For transverse propagation with the electric 
field perpendicular to the dc magnetic field (E^=0) , one 
obtains ^ ^ 

.2- ,2 V- 

- “’"o -k 7~ = -'"o — 


1 


1 


( 2 - 22 ) 
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This node is generally elliptically polarized in a 
plane perpendicular to the .do magnetic field. It 
should be noted that the propagation constant for the 

r 

transverse direction is reciprocal, that is, reversal 
of the direction of the propagation (or alternatively, 
the direction of the dc magnetic field) does not 
change the propagation constant. However, the 
reflection coefficient for a wave incident on the 
m.agnetoplasma where the electric field is in the plane 
of incidence is non-reciprocal, and this is considered 
for the development of non-reciprocal devices for use 
at millimeter and submillimeter wavelengths in Part II. 



PART I 


NON-RECIPROCAL WAVEGUIDE ISOLATORS 
FOR USE AT MILLIMETER WAVELENGTHS 

The purpose of Part I is to present a theoretical 
and experimental study of the design of a new semi- 
conductor waveguide isolator at 94 GHz where the 
atmosphere has a window and attenuations are frequently 
as low as 1 or 2 db . In Chapter 3 a theoretical 
analysis of propagation of em waves in a circular wave- 
guide containing an annular column of a solid state 
magnetoplasma is given. First a heuristic argument 
based on the exclusion or absorption of the em waves 
due to the magnetoplasma is presented. Second a mode 
matching analysis (henceforth abbreviated M4A) is 
employed in the theoretical analysis. Chapter 4 
presents the experimental apparatus and procedures 
used to measure the characteristics of propagation in 
the circular waveguide containing the annular column 
of the solid state magnetoplasma. The properties of 
the material used in the experiments are also given. 

In Chapter 5 the experimental results are given along 
with the theoretical results from the approximate 
theory and from MMA. The correlation between the 
expe;rimental and theoretical results and the validity 
of the assumptions made in the approximate theory are 
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discussed. Chapter 6 presents the sutninary of the work 
of Part I and conclusions. 



Chapter 3 


PROPAGATION CHARACTERISTICS OF CIRCULAR WAVEGUIDE 
CONTAINING AN ANNULAR t^GNETOPLASMA COLUMN 

The propagation characteristics of em waves in a 
circular v/aveguide- partially ’ filled with a plasma have 
been investigated extensively examining the effect of 
a plasma upon waveguide m.odes . For example, R.N. 

Carlile [1966] measured the propagation constants of 
the TE^^ modes in a cylindrical waveguide partially 
filled with a gaseous plasma in a longitudinal dc 
magnetic field and observed that right and left cp 
waves travel with different propagation constants. 

The theoretical and experimental study of a circular 
35 GHz waveguide containing a coaxial InSb rod was 
done extensively by B.R. McLeod and W.G. May [1968, 

1971] . 

The purpose of this chapter is to' present a 
general theoretical analysis of the characteristics 
of microwave propagation in a circular waveguide 
containing a finite length of an annular plasma column 
with a longitudinal dc magnetic field, and a more 
specialized theoretical study of a 94 GHz waveguide, 
isolator employing this geometry. The basic operation 
of the isolator depends on the exclusion or absorption 
of cp waves by the plasma. The extent of the exclusion 
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or absorption depends on the sense of the circular 
polarization of the eiti wave with respect to the vector 
direction of the longitudinal dc magnetic field. 


3~'l Effect of Nearly Zero Permittivity on Circularly 

P olarized Waves 

With a longitudinal dc magnetic field on the 
plasma column, the total dielectric constant of the 
material due to both lattice and free carriers seen by 
the em v;aves depends on the direction of the circular 
polarization with respect to the longitudinal magnetic 
field. The total relative dielectric constant derived 
in Section 2-2 has the form: 



(3~1) 

(3-2) 

(3-3) 


The subscript ± refers to the right and left cp 
components in the direction of the dc magnetic field. 

It is readily seen that the real part of the dielectric 
constant could be relatively large for a microwave 
signal circularly polarized one way but very low or 
even negative for the opposite cp wave and the imaginary 
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part small in both cases. . This characteristic was 
utilized for the development of the non-reciprocal 
devices at millimeter wavelengths by exciting a right 
or left cp em wave in the circular waveguide containing 
the annular plasma column as discussed in the next 
section . 

3-1-1 Approximate Theory for Attenuation 
Coefficient 

Consider the configuration shown in Fig. [3-1]. 
Using a heuristic argument let us estimate the ratio 
of the electric field strength in the annular plasma 
column to the field in the center section of air- 
filled region. It is assum.ed that the fields both in 
the center section of the air-filled region and in the 
annular plasma column are of mode pattern. The 

amplitude of field in each region can be, however, 
drastically different and depends on the dielectric 
constant of each region. Of course, the above 
assumption is not truly valid and will not be used in 
MMA in the next section. 

Using a simple static analysis the ratio of the 
field concentrations is 



(3-4) 
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where K is the relative dielectric constant of the 

solid state magnetoplasma, and E. and E are the 

X o 

uniform fields in the air-filled center section and in 
the annular column, respectively. 

If the dc magnetic field is adjusted so that as 

in Eq. (3-1) is close to zero, it appears reasonable 
from Eq. (3-4) that the left cp waves are strongly 
excited in the annular plasma column. Thus, the wave 
propagates almost exclusively inside the annular 
plasma column, and because of the lossy part of the 
dielectric constant of the plasma, much of the left 
cp wave will be absorbed by the plasma and dissipated 
as heat. On the other hand, the same em wave traveling 
from the opposite direction v^ith respect to a 
longitudinal dc magnetic field is right circularly 
polarized. It is predicted from. Eq. (3-4) that the 
right cp wave is largely excluded from the annular 
plasma column because of large (~32) and is only 
slightly absorbed by the plasma. 

Let us examine the attenuation coefficient using 
the approximate model discussed above. Consider the 
infinitely long annular plasma column placed in a 
uniform em field with the propagation direction 
parallel to the axis of the column. Assuming that the 
transverse fields in the annular plasma column are 
approximately uniform for both polarizations -of the em 



waves f the power dissipation per unit length in the 
plasma annulus is 
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^diss 



(3-5) 


where E is the field inside the annular plasma column 
and is the area of its cross section. From the 
foregoing discussion, the field inside the annular 
plasma column is also assumed to be a uniform TE mode 
and the power flow inside the annular plasma column 
is then approximately given by 


P 

o 


1 

2 


' o ' , _ 1 e ± 

ri o 2 Z 
o o 



] , 

(3-6) 


where is the impedance of the TE mode in the annular 
plasma column. Sim.ilarly the power flow in the center 
section of the air-filled region P^^ is given by ■ 


P. 

1 




(3-7) 


o 

where E^ is the field in the air-filled center section 
and is the area of its cross section. By neglecting 
the reflections at the plasm.a boundary, the approximate 
attenuation coefficient a for both senses of the 
circular polarization of the era vraves due to absorption 
in the plasma is found (by dividing the dissipated 
power by the total power flow) to be 



1 e ± 


1 we K"|e j^A 

2 o ± ' o ‘ o 





o 


• (3-8) 
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Using Eq. (3-4) to eliminate the and one obtains 

for the approximate attenuation coefficient: 


a . = 


we 2 K" |1+K i^A 
o O ± ' ± ' O 


2 . 

ji 


R. 1 1+K^ 1 2 [1- (‘^) ] a^+4 I 

c 


(3-9) 


e ' ‘ 


'A. 


This expression is used to estimate the attenuation co- 
efficient in the annular plasma column and is compared 
with that from MIA discussed in the next section. 


3-2 M ode Matching Analysis to Boundary Value Problem 

In the previous section, it is assumed that the 
field configuration is of mode pattern even in the 

presence of the plasma. However, such a simplification 
is not truly applicable and will not be used in the 
present analysis, because a single incident mode can 
excite infinite number of gyromagnetic modes at the 
ends of plasma column. These modes in turn excite an 
infinite number of TE and TM modes in the empty 
circular waveguide. Most of these modes are evanescent 
but even so must be included along v/ith the propagating 
modes- in the rigorous treatment of the boundary 
conditions . 

The general calculations of the reflection and 
transmission coefficients of a finite-length section of 
the annular magnetoplasma column consist of the follow- 
ing two steps ; 
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1) Determination of the propagation constants of 
the gyromagnetic modes in the section of the 
annular magnetoplasma column, OiziS-. 

2) Use of these propagation constants to match 
the boundary conditions at the ends of 
plasma column , z—Oft. Since the gyromagnetic 
modes in the section of the annular magneto- 
plasma column do not form orthogonal sets^ 
the point-matching method is employed to 
determine the reflection and transmission 
coefficients at each boundary, z=0,£- 

These steps will be performed in Section {3-2-1) and 
(3-2-2) , 

3-2-1 Characteristic Equation for Waveguide with 
Annular Homogeneous Magnetoplasm.a -Column 

The waveguide system considered here is illustrated 
in Pig. [3-1]. It consists of an annular homogeneous 
plasma column with axis parallel to a dc magnetic field 
in the z-direction. 

In terms of the complex tensor 'dielectric constant 
K given in Eq. (2-5) and for the fields that vary as 
expCjmt)^ Maxwell ' s equations inside the annular plasma 
coliomn are^ given by 

VxE = 


(3-10) 
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VxH = iwe K«'E , 
o 

and 

V* (e t-E) = 0 f 
o 

V- (y^H) = 0 . 


{3"1D 


(3-12) 

(3-13) 


Since the solutions of the eiti fields satisfying the 
Max\irell’s equations have a z-dependence of the form 
exp(-rz), where the propagation constant F is indepen- 
dent of Zf it is convenient to make following separa- 
tions : 


E 



a E 
z z 


(3-14) 


H = Hrr, + a B 
T z z 


(3-15) 


V = 


T 


+ 


* 3 

^z 3z 



S,(-D . 


(3-16) 


The coupled wave equations describing the longitudinal 
components ,of the em waves in the region of the annular 
homogeneous plasma column aSrSb are given by A. Bers 
[1963] ; 


vJe + aE 
T z z 

= bH , 
z 

(3-17) 

V^H + cH 
T z z 

= dE , 
z 

(3-18) 


where 


= (F^+k^Kj^) 


"1 


b = 


jwy^r 


5 

K 


1 
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r2*v2 V- 

-rr 




d = -jwe^r 




K 


1 


This coupled set of second-order equations can be trans- 
formed into an uncoupled set of fourth-order equations : 

[V^ + {a+c)V^ + (ac-bd)]E^ = 0 , (3-19) 

and 

[V^ + (a+c)vj + (ac-bd)]H^ = 0 . (3-20) 

By admitting solutions of the form exp(-jp’r^) for the 
transverse components of the em waves in the region of 
the annular plasma column a<rSb, the dispersion relation 
is given by 

p^ - (a+c)p^ 4- (ac-bd) = 0 . (3-21) 

The longitudinal components of the electric field can 
be obtained from the following differential equation: 

(Vj + P?)E^. = 0 . (3-22) 

The acceptable solution which describes the longitudinal 
component of the electric field inside the plasma 
(airib) is expressed by the linear combination of 
ordinary Bessel functions and Neumann functions : 

(3-23) 

The longitudinal com.ponent of magnetic field is given 
by 
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+h2{CJ^{p2r)+DN^(p2r) }] 


( 3 - 24 ) 


where 


h 


1,2 


2 

^"Pl,2 


c-p 


1,2 


The transverse components of the electric and magnetic 
fields can be derived after lengthy tedious computations 
and are given in Appendix B. 

The wave equations in the center section of the 
air-filled region (r<a) are 


and 


where 


= 0 


V^H° + = 0 


2 p2 ^ , 2 

a = r + k 

o 


( 3 - 25 ) 


( 3 - 26 ) 


Two independent solutions which describe the longitudinal 
components of the electric and magnetic fields in the 
center section of the air-filled region are 


E = 


DJ (qr)e 
m ^ ' 


jm(})-rz 


( 3 - 27 ) 


and 


H° = VJ 
z m ^ 


( 3 - 28 ) 
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Using 


F° ~ 


->-o 

= 


q q 


r n ^ -t gO 

2Vz 2 z Vz ' 

q q 


(3-29) 

(3-30) 


all transverse components of the electric and magnetic 
fields in the center section of the air-filled region 
are listed in Appendix B . 

The boundary conditions along a radius are that the 
tangential components of the electric fields are zero at 
r=b and that the tangential components of the electric 
and magnetic fields are continuous at r=a. The character- 
istic equation turns out to be a 6x6 determinant and is 
given in Appendix B. 


3-2-2 Matching Fields at the Ends of the 
Magnetoplasma Column 

For a plasma section of finite length, one is con- 
fronted with the problem of matching fields at the ends 
of the plasma column. This step must be taken in order 
to derive matrix reflection and transmission coefficients 
at the ends of plasma column. The over-all attenuation 
due to the finite section of the annular plasma column 
is then calculated from the transmission coefficient of 
the dominant mode at the end of the plasma column as 
given later in Eq. (3-43) . 
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The problem of matching fields at the ends of the 
plasma column is quite complicated for the following two 
reasons. First, since a single incident mode can excite 
a number of gyromagnetic modes in the plasma section and 
these, in turn, excite a number of TE and TM modes in 
the empty waveguide, these higher modes must be included 
in applying the boundary conditions . Second , since the 
gyromagnetic modes of the plasma section do not form 
orthogonal sets due to the existence of loss in the 
plasma [Clarricoats , 1970], the point matching method 
is introduced to compute the reflection and transmission 
coefficients at the ends of the plasma column. 

Considering the configuration shown in Fig. [3-2] , 
the continuity of the tangential electric and magnetic 
fields at z~0 gives 


and 










(3-31) 

(.3-32) 


Similarly the continuity of tangential electric and 
magnetic fields at z=a leads to 


[A|°h[E^^°h + [^°^] = [a|^^][E 


T(l) 




f 


(3-33) 
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and 

+ [b|^^3 . (3-34) 

Here [E] and [H] are column matrices of electric and 
magnetic fields, and the superscripts T and R indicate 
transmitted and reflected em waves, respectively. The 
index of superscript (0) indicates the field in the 
empty waveguide and (1) indicates the field in the 
plasma section. The subscript t of E and H stands for 
the transverse component of em waves. [A] and [B] are 
the row matrices of the amplitude coefficients of the 
transmitted and reflected em fields, respectively. The 
subscript 0 of A and B indicates the transmission and 
reflection coefficients at z=0 and the siabscript I 
indicates the transmission and reflection coefficients 
at z=a. The explicit forms of the transverse components 
of the em fields in the magnetoplasma section and in the 
empty waveguides are given in Appendix C. 

For the case when a single dominant mode is incident 
at 2=0 and no reflections take place at z=®°, 

rA^°^3 = (1,0,.. .,0) , (3-35) 

and 

= ( 0 , 0 ,. .., 0 ) . 


■ (3-36) 
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The row matrices of other reflection and transmission 
coefficients are given by 


ra(0), _ /rp(h,0) T,(h,0) „{e,0] rn(e,0)v 

LA^- i - **-^nt ' It ** lat ’ 


( 3 - 37 ) 


( 0 ), 

o 


, (h,0) j.{h,0) j.(e,0) R(e,0). 

^^It *’*^nt ' ^It * * mt ^ 


( 3 - 38 ) 


lA^^h = (T^^^ 




^(1) ) 
'^n+m,t^ " 


™(i) ~r . 


( 3 - 39 ) 

( 3 - 40 ) 


[e'1>1 = (R<1> 


) 

n+m,t^ 




( 3 - 41 ) 

( 3 - 42 ) 


Here T and R indicate the transmission and the reflection 
coefficients. The superscripts (h^O) and (S;.0) indicate 
the TE mode and the TM mode in the empty waveguide, 
respectively, and (1) indicates the gyromagnetic mode 
in the magnetoplasma section. The subscripts l....m,n 
indicate mode numbers , and t stands for the transverse 

component of em waves, is the propagation constant 

‘t'll ‘ 

of the m mode in the magnetoplasma sectxon. 

These reflection and transmission coefficients are 
evaluated by numerical techniques using the point 
matching method. The boundary conditions are imposed 
at only a finite nxiraber of points along the contour of 
the cross section. This formulation becomes increasingly 
exact by increasing the number of modes considered and 
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the number of points chosen. 

The over-all attenuation due to the finite section 
of the annular plasma column is finally given by 

Loss = -20 log db (3-43) 



Chapter 4 


EXPERIMENTAL PROCEDURE 

This chapter presents the experimental apparatus 
and procedures in the experimental investigations for 
the circular waveguide containing the annular plasma 
column in a longitudinal dc magnetic field. The 
characteristics of the semiconductor used as a plasma 
medium will be established. Discussions are extended 
to the sample preparations and the experimental proced- 
ure in observing the non-reciprocal phenomena in the 
plasma waveguide. 

4-1 Properties of Solid State Plasma 

The material chosen for the solid state plasma is 
n-type InSb, which has a high mobility of free carriers 
and a relatively low carrier concentration. The average 
value of the effective mass for n-type InSb at 77 °K from 
infrared Faraday effect measurements was m*=0.0143m^ 
[Hogarth^ 1965] . From the Faraday rotation experiments 
on ,n-type InSb at 77 °K performed in this laboratory by 
B.R. McLeod [1968], the static relative dielectric 
constant was found to be 16.0. The mobility of the 
free carriers and the carrier concentration are deter- 
mined by performing either Hall measurements or 
van der Pauw measurements [van der Pauw, 1958] . For 
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nondegenerate semiconductors and small magnetic fields 
the Hall coefficient is given by 



r 2 
^ (P+bn) ^ 




where 



For n>>p, the Hall coefficient reduces to 



(4-1) 


(4-2) 


The factor r^^ is determined by the scattering mechanism 
[Smith, 1964] . _ When acoustic scattering is dominant, 
r =1,18. If ionized impurity scattering dominates, 
rj^=1.93. However, the scattering of conduction 
electrons by both impurity ions and the lattice are 
dominant at 77 °K in n-type InSb. V.A. Johnson and 
W. J. Whitesell [1953] , and R.T. Bate, R.K. Willardson 
and A.C. Beer [1959] have made an extensive theoretical 
and experimental study of magnetoresistance effects in 
type InSb between 50 °K and 200°K, and the value r^^ 
was found to be 1.09. After making this correction, the 
parameters of n-type InSb used in the experiment were 
found as listed in Table [4-1] . 
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4-2 Experimental Set-up 

The devices designed for non-reciprocal phenomena 
consist of two transitions, two polarizers and the test 
section shown in Fig. [4-1] . TE^q mode propagating in 
the rectangular waveguide excites a linearly polarized 
TEii mode in the circular waveguide due to the transi- 
tion. The polarizers which consist of quarter-wave 
teflon slabs in the circular waveguide set at 45° with 
respect to the electric field in the rectangular wave- 
guide convert the linearly polarized TE^^ mode into the 
cp TE^j^ mode. The test section contains the annular 
magnetoplasma column of n-type InSb. These columns of 
inSb were cut by using an ultra-sonic drill and polished 
chemically in a cromine and methanol etch. They were 
supported-by styrofoam in the test section. The small 
opening between the column and the circular waveguide 
wall was filled with the silver paint. 

The devices were placed in a longitudinal dc 
magnetic field using a 6-inch Varian magnet with a 
2-inch gap. A styrofoam Dewar surrounding the semi- 
conductor was filled with the liquid nitrogen. The 
microwave system was composed of a klystron, an 
isolator, a variable attenuator and a crystal detector. 
The klystron provided a 94 GHz signal modulated at 
1 kc/s. A lock-in amplifier was used to take the data. 
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A block diagram of the experimental set-up is shown in 
Fig . [4-2] . 


TABLE 4-1 

Electrical Properties of n-type InSb Crystal 


Electron 
Effective Mass 

Dielectric 

Constant 

Electron 

Mobility 

Carrier 

Density 


75°K 

0.0143m 

o 

16.0 

5.2xl0^cm^/v-sec 

2.8xl0^'"/cm2 




Fig. 4-1. Waveguide Isolator containing Annular Plasma Column. 




Fig. 4-2. Block Diagram of Experiment Set-up. 











Chapter 5 


DISCUSSION OF EXPERIMENTAL AND THEORETICAL RESULTS 

This chapter presents the experimental and theoreti- 
cal results of propagation characteristics for a circular 
waveguide containing an annular InSb column in a longi- 
tudinal dc magnetic field. The results calculated from 
the approximate theory of Section 3-1 and from the 
MMA of Section 3-2 are given. The validity of the 
assumptions made and the correlation between the 
experimental and theoretical results are discussed. 

5-1 Experimental Results 

Fig. [5-1], [5-2] and [5-3] show the experimentally 
determined attenuation for the circular waveguides 
containing the various lengths I of InSb. Only losses 
within the test section are included. Attenuation could 
be measured to within ±0.5db and the magnetic field was 
measured with a gaussmeter to an accuracy of ±100 gauss. 

5-2 Comparison with the Approximate Theory 

The calculated attenuation curves derived from the 
approximate theory of Section 3—1 are shown in Fig- 
[5-4] . This figure shows that approximate theory gives 
a fairly good agreement with the corresponding 
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Fig. 5-1. Experimental Results for Hollow InSb Waveguide 
Isolator [Parameters of InSb are given in 
Table 4-1. Geometry of 94 GHz Waveguide 
Isolator is shown in Fig. 4-1 with 
^.=0 . 127mm] . 
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Fig. 5-2. Experimental Results for Hollow InSb Waveguide 
Isolator [Parameters of InSb are given in 
Table 4-1. Geometry of 94 GHz Waveguide 
Isolator is shown in Fig. 4-1 with 
£=0 . 762mm] . 
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Fi.g. 5-4. Theoretical Results for Hollow InSb Waveguide 
Isolator (5. — 0.7 62mm) Based on Approximate 
Theory [Parameters of InSb are given in 
Table 4-1. Geometry of 94 GHz Waveguide 
Isolator is shown in Fig. 4-1 with A=0. 762mm] 
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experimental results for a low longitudinal dc magnetic 
field up to approximately 6 kG. The assumotion 
of a uniform plane wave propagating in the section of 
the annular InSb column appears to be quite reasonable 
and the exclusion or the absorption of em waves is the 
physical explanation summarized as follows - Since the 
power density of the right cp v/ave is very small in the 
InSb and is heavily concentrated in the air-filled 
center section, its propagation is virtually independent 
of the longitudinal dc magnetic field. On the other 
hand, since the left cp wave is excited stronaly i-^ithin the 
lossy InSb, the em waves are absorbed and dissipated 
as heat. For the case of higher longitudinal dc 
magnetic fields, MMA of the boundary value problem 
shows that the excitation of higher gyromagnetic modes 
becomes significant, and some deviation of the 
electric field from the mode pattern is expected. 

Thus this causes some trouble in the approximate theory. 

It is also found that the reflection coefficient for 
the dominant mode obtained from MMA is relatively small 
(<0.5) at the end of plasma column, and thus it is a 
reasonable first approximation to neglect the attenua- 
tion due to the reflection. But the approximate theory 
did not account for the transition between the empty 
waveguide and the guide in which the annular InSb 
column was mounted, nor did it account for my reflections, 
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and thus extremely good agreement between theoretical 
and experimental results was not expected. 

5-3 Comparison with Mode Matching Analysis 

The first step of MMA was to determine the propaga- 
tion constants of the gyromagnetic modes. Since the 
characterisjti'C equation is quite complicated, an intera- 
tive root-finding technique was employed to find the 
propagation constants. The real and imaginary parts of 
the propagation constants for both right and left 
polarizations were obtained by computer and are shown 
in Fig. [5-5-] , [5-6] . Three place accuracy was used 
throughout the 'calculations . 

The resulting propagation constants are utilized 
to match the boundary conditions at the ends of the 
plasma column. This second step deteptiines the 
amplitudes ' of the reflection and transmission coeffic- 

^ I 

ients of the. gyromagnetic modes and of the empty wave- 
guide modes. In the present case, four gyromagnetic 
modes in the section of the annular InSb column and two 

. s' 

TE and two TM modes in the empty waveguide are taken 
into account. Point matching was employed to solve the 
inhomogeneous linear equations of order, 16 . The overall 
attenuation due to the various lengths of the InSb 
columns was evaluated from' the transmission coefficients 
of circular mode at the end of the column and is 
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Fig. 5-5. Noriralized Propagation Constant for Right Circularly Polarized 
Wave as a Function of Magnetic Field [Parameters of InSb are 
given in Table 4-1; Geometry of 94 GHz Waveguide Isolator is 
shown in Fig. 4-1] . 
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Fig. 5-6. Normalized Propagation Constant for Left Circularly Polarized 
Wave as a Function of Magnetic Field [Parameters of InSb are 
given in Table 4-1. Geometry of 94 -GHz Waveguide Isolator 
is shown in Fig. 4-1] . 
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shown in Fig. [5-7 ] , [5-8] and [5-9] . The general 
agreements between the eKperimental and the theoretical 
results are quite good. The boundary conditions are 
imposed at only four points at the ends of the plasma 
column by taking into account only four gyromagnetic 
modes in the magnetoplasma section and four higher 
excited modes in the empty waveguide section. It is, 
however, expected that, although the inclusion of higher 
excited modes as well as the dominant modes is essential 
in a rigorous treatment of the boundary conditions at 
the ends of the plasma column, the convergence is fairly 
rapid [Champlin, et al., 1969]. Thus small discrepancy 
between theoretical and experimental results appears to 
be mainly due to experimental error as discussed as 
follows . The annular InSb column was cut by using an 
ultra- sonic drill. Although the resulting InSb column 
was close to an annular cross section, there was some 
departure from the circular symmetry required for a 
circularly polarized wave, and this caused some trouble 
in the experiments. The results of the van der Pauw 
measurements showed variations by as much as 10% 
existed in the crystal properties. This inhomogeneity 
would also explain some of the discrepancy between 
theoretical and experimental results . 

The reflection coefficients of the incident TE^^^ 
mode at the ends of the plasma column show that the 
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Fig. 5-7. Theoretical Results for Hollow InSb Waveguide 
Isolator [Parameters of InSb are given in 
Table 4-1. Geometry of 94 GHz InSb Waveguide 
Isolator is shown in Fig. 4-1 with ^=0,1273ran]. 
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Fig. 5-8. - Theoretical Results for Hollow InSb Waveguide 
Isolator [parameters of InSb are given 
in Table 4-1. Geometry of 94 GHz InSb 
Waveguide Isolator is shown in Fig. 4-1 
with £=0. 762mm]. 
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Fig. 5-9.. Theoretical Results for Hollow InSb Waveguide 
Isolator [Parameters of InSb are given 
in Table 4-1. Geometry of 94 GHz InSb 
Waveguide Isolator is shown in Fig. 4-1 
with £=1. 524mm]. 






55 


devices are absorbing the, microwave power to achieve 
the attenuation obtained experimentally. A small 
amount of the reflected wave is linearly polarized in a 
horizontal plane/ and thus is absorbed by the resistance 
card placed behind the polarizer. Of course, it could 
be possible to prevent this reflection by a careful 
design of customary matching transformers such as a 
conical taper or an end cap of dielectric material 
whose dielectric constant is chosen so as to provide an 
impedance match at the ends of plasma column. 



Chapter 6 


SUJffilARy AND CONCLUSIONS OF PART I 


Theoretical and experimental investigations were 
performed for the purpose of the utilization of non- 
reciprocal behavior in the solid state magnetoplasma 
for the development of millimeter waveguide isolators. 

The geometry employed consists of a circular waveguide 
with a coaxial annular column of a semiconductor plasma 
in a longitudinal dc magnetic field. 

The basic operation of this type of isolators is 
due to exclusion or absorption of em waves in a magneto- 
plasma. The microwave signal for the right cp wave is 
excluded from the magnetoplasma, and thus the signal can 
be only-slightly absorbed. The microwave signal for the 
left cp wave, however, penetrates into the lossy magneto- 
plasma and is strongly absorbed. The approximate theory 
is based on this physical reasoning. A mode matching 
analysis is based on the solutions of the boundary 
value problem- The propagation constants of the gyro- 
magnetic modes were determined in the section of a 
circular waveguide with an annular column of magneto- 
plasma. For the, case of a waveguide with finite length 
of an annular plasma column, the boundary conditions 
were satisfied at the ends of plasma column by taking 
into account the higher excited modes as well as the 
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dominant mode. The reflection and transmission 
coefficients were determined at the ends of the plasma 
column by using the point-matching method. 

The experiments were performed at 94 GHz. In 
order to obtain characteristics useful for an isolator, 
i.e. low forward loss and high reverse loss, the InSb 
isolator was operated at 75 °K. In those systems in 
which cooling is already present such as in maser 
systems, this is not a serious problem. The experiments 
demonstrated 23 db isolation with 7 db forward loss for 
one waveguide isolator and corresponding theoretical 
data calculated from a mode matching analysis were 
36 db isolation with 5 db forward loss. - This dis- 
crepancy between theoretical and experimental results 
mainly results from the incomplete theoretical analysis. 
It is assumed in the approximate theory that the field 
both in the center section of air-filled region and in 
the annular column are of TE^^ mode pattern, and the 
amplitudes of the electric fields are determined by the 
dielectric constant of each section. However, it is 
found from a mode matching analysis that some deviation 
of electric field from the TEj^^ mode pattern in these 
regions takes place, and this causes some trouble in 
the approximate theory. In the mode matching analysis 
the point matching method was adopted to solve the 
boundary problem. The boundary conditions are imposed 
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at only four points at the ends of plasma colmn by 
taking into account four gyromagnetic modes in the 
magnetoplasma section as well as four excited higher 
modes in the empty waveguide region. It is possible, 
of course, to correct these by increasing the number of 
modes considered and the number of points chosen. 

The hollow InSb waveguide isolator can be operated 
with a wider band width, since the maximum attenuation 
was insensitive to the dc magnetic fields (see Fig. 

[5-1] , [5-2] and [5-3]) and hence to the operating 
frequencies . It is speculated that the lower forward 
loss could be obtained in an annular InSb column with a 
larger diameter hole but with expense of isolation. 

Although the requirement of cooling may limit the 
applications for millimeter wavelength devices, the 
solid state magnetoplasma finds important applications 
for the development of the non-reciprocal devices in 
the submillimeter wave region. 



PART II 


NON-P.ECIPROCAL REFLECTION BEAI'l ISOLATOR 


FOR USE AT MILLIMETER AND 
SUBMILLII4ETER WAVELENGTHS 


The purpose of Part II is to, discuss the non- 
reciprocal reflection from solid state magnetoplasmas 
for the development of the reflection beam isolator for 
use at millimeter and submillimeter wavelengths. 

Chapter 7 presents the theoretical deviation of reflec- 
tion coefficients’ from a plasma medium. In Chapter 8 
the experimental apparatus and procedures used to 
measure the reflection coefficients from the plasma 
interfaces are given. The properties of the semi- 
conductors used as plasma media are also given. The 
brief discussion on the design and construction of HCN 
laser used as a far-infrared source for the experiments 
and on the IR detectors are also given. Chapter 9 
presents the theoretical and experimental results . The 
results for the reflection from the interface between a 
dielectric medium and GaAs at 94 GHz, and from the inter- 
face between free space and dielectric coated InSb at 
the wavelength of 337 ym are presented. In Chapter 10 
physical explanation of the non-reciprocal reflections 
is presented. Possible refinements for the applications 
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of the devices and comparison with other devices are 
also brought out. 



Chapter 7 

THEQRETICAL DERIVATION OF REFLECTION COEFFICIENTS 


Many studies have been made on the nature of the 
reflection from the ionosphere. N.F. Barber and D.D. 
Crombie [1959] have calculated the reflection coeffic- 
ient from a sharply bounded ionosphere in the presence 
of a horizontal magnetic field perpendicular to the 
plane of incidence and found that the reflection 
coefficient for em waves in the ionosphere incident 
from the west is' niamerically greater than that for the 
wave incident from the east. Non-reciprocal reflection 
of em waves incident on a solid state raagnetoplasma at 
94 GHz was studied by J.M. Seaman [1968] . 

This chapter outlines the theoretical analysis of 
the reflection coefficient of a wave incident on a 
magnetoplasma, which depends on the orientation of the 
dc magnetic field and the direction of propagation. We 
consider only the case in which the plane of polariza- 
tion of em waves is in the plane of incidence , and the 
direction of propagation is perpendicular to a dc 
magnetic field which is parallel to the surface of the 
magnetoplasma (Kerr transverse magneto-optical effect) . 

Three configurations are considered with the goal 
of achieving best isolation with low insertion loss. 
These are the reflections of a plane wave from the 
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interfaces between free space and a magne^oplasma, between 

I ^ * ' ' » ' ' t ■'if « I \ 

a dielectric medium and a magnetoplasma, and between 
free space and a dielectric coated magnetoplasma, as 
shown in Fig.' [7-1]. 

Define Cartesian coordinates where a dc magnetic 
field is applied in the z-direction. For a plane wave 
with harmonic time dependence, exp (jut). Maxwell's 
equations are given by 


-jwy^H = VxE , 

(7-1) 

jus E = R"^VxH . 

(7-2) 


1“^ is the inverse of the relative dielectric tensor, 
of a plasma and is given in the following form: 

A -D 0 

d' a 0 (7-3)' 

0 0 B 
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1-1 Reflection from Air-Plasma Interface 

The reflection coefficient R for a plane boundary 
between free space and a plasma medium (see Fig. [7-la]) 
has been calculated [Wait, 1962] and is given by 


COS0 - A 
COS0 + A 


where 


A = A‘ 


(i - 


sin^0) 


+ D • sin0 


(7-5) 


(7-6) 


7 “2 Reflection from Dielectric-Plasma Interface 

The calculation is extended below to include the 
effects of a dielectric medium placed on the plasma as 
shown in Fig. [7-lb] , and the reflection coefficient R^^ 


from the interface between a dielectric medium and a 
plasma medium is given by 

cose - A. 


_ M 

^ ~ COS0 + A,, ' 

M 


(7-7) 


where 


3s 


and 




A^ = Kj^EA- (^ - sin^0) + D-sin0] , (7-8) 

is a relative dielectric constant of a dielectric 


medium. 
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7-3 Reflection from Air-Dielectric Coated Plasma 

Interface 

In this section the reflection coefficients from 
the. interface between free space and a dielectric 
coated magnetoplasma are derived by both the field 
matching method and the impedance method. 

7 - 3-1 Field Matching Method 

Generally the magnitude of all field components of 
em waves reflected from a multiple layer interface can 
be evaluated by matching the tangential components of 
the field at each plane interface. 

Considering the configuration as shown in Fig. 

[7-2], where the electric ‘field is polarized in the plane 
of incidence and the fields do not vary in the z 
direction, we are required to solve a 4 order inhomo- 
geneous linear equation in order to satisfy the boundary 
conditions. Although it is quite tedious to solve the 
inhomogeneous linear equation, the method is especially 
useful to evaluate magnitudes of all components of em 
waves inside a plasma medium. These components are used 
in Chapter 10 to calculate the ellipticity of the 
electric field and the orbit of electrons in the magneto- 
plasma in order to study the physics behind the non- 
reciprocal reflection. 
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Fig. 7-2. Orientations of Fie^d Vector E, 
Propagation Vec;^or V and dc 
Magnetic Field B. 
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The boundary conditions are that the tangential 


components of 

the 

fields 

are 

; to be continuous 

at 

y=0. 

S, in Fig. [7- 

2] . 
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(7-12) 


Here E and H indicate the electric and magnetic fields, 
and the corresponding superscripts t and r stand for 
transmitted and reflected em waves. The right sub- 
scripts X and z indicate the x and z components of 
fields f and’ the left s\ibscripts 0, 1 and 2 indicate the 
fields in free space, in the dielectric and in the mag- 
netoplasma. T and R indicate the transmission and ref- 
lection coefficients, and the corresponding subscripts 
0, 1 and 2 indicate amplitude coefficients in free space, 
in the dielectric and in the plasma. The detailed cal- 
culations of all field components in free space, in a 
dielectric, and in a magnetoplasma are given in the 
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standard textbooks [Collins,, 1960; Stratton, 1941]. 

The reflection coefficient from the interface 
between free space and a dielectric coated magneto- 
plasma, and the amplitude of the field components in 

t ti b 

the magnetoplasma, „E and , are obtained by 

. 4(!!# Z ^ ^ j 

'til • 

solving the 4 order inhomogeneous linear equation as 
given in Eg. (7-9) to (7-12). 


7.3,2 Impedance Method 


The wave impedance G for the transmitted wave is 


defined by 


,t 


. E 
H 


t 

X 

t 

2 


(7-13) 


and for the reflected wave, 

T 

E 

Vj ^ ; 


(7-14) 


z 

where superscripts t and r indicate the wave impedances 
for transmitted and reflected waves, respectively. The 
wave impedances G*s for each layer shown, in Fig. [7-2] 
are as follows: In free space (0) : 


g'*^' = = Z COS0 . (7-15) 

0 o o • 

In a dielectric medium (1)*; 

4- 2 h 

G^ = Gf = ~ (K„ - sin^e)- V (7-16) 

1 1 K,, M 

M 

where Kj^ is the relative dielectric constant of the 
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dielectric medium. In a magnetoplasma (2) : 

h 

gJ = Z [A* (J - sin^e) + D-sin6] (7-17) 

h 

= Z [A- - sin^e) - D-sin0] (7-18) 


Here 0 is the angle of incidence (measured to the 
negative y axis in free space as shown in Fig. [7-2]). 
Note that because of the quantity D due to an anisotropic 
medium, G 2 and G 2 are not equal as they would be in 
the isotropic medium. This non-reciprocal wave imped- 
ance is explained physically by the excitation of em 
waves at the dielectric-magnetoplasma interface due to 
the oblique incident em waves, which is given in 
Section ‘10-1 in detail- This is the one of the key 
factors for the explanation of the non-reciprocal 
reflection. 

The reflection coefficients at the interface 

t 

between a dielectric medium and a plasma medium 
in Fig. [7-2]) is, therefore, given by 


and 



(7-19) 


(7-20) 
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From analogy with transmission line theory^ the 
impedance seen at the interface between free space 
and a dielectric medium (y=0 in Fig. [7-2]) is given by 


Z 


1 



1 + r^^exp (-j2k^£,j^) 

1 + R^exp {- j 2kj^5.^) ' 


(7-21) 


2 

where k^ = k^ (K^^ - sin 9 ) is a wave number in a dielec- 
tric medium and 5-^ is the thickness of the dielectric 
layer on the plasma medium. 

Finally the reflection coefficient at the interface 
between free space and a dielectric coated magnetoplasma 
(y-0 in Fig. [7-2]) is given by 



(7-22) 


If D given in Eq. (7-17) and (7-18) is complex due to 

the presence of collisions in a magnetoplasma, it will 

be found that the reflection coefficient R is non- 

o 

reciprocal when a dc magnetic field is reversed. That 
is,' reversal of the direction of propagation, dr alter- 
natively the sense of the dc magnetic field, changes the 
magnitude of the reflection coefficient. The physical 
explanation for the necessity of collisions in the non- 
reciprocal reflection is given in Section 10-1 in 
detail. This phenomenon is exploited to develop non- 
reciprocal devices such as isolators for use at milli- 
meter and submillimeter wavelengths. 



Chapter 8 


EXPERIMENTAL PROCEDURE 


This chapter is concerned with establishing the 
electric properties of the semiconductors used as the 
plasma medi'um. The experimental apparatus and the 
procedure for the measurement of the reflection from 
the magnetoplasma are presented. Brief discussions 
on the design, construction and operation of HCN laser 
and' IR detector are also given in this chapter. Details 

t 

are presented elsewhere [Kanda, 1971] . 

87 I Properties of Solid State Plasma 

At first the experiments were demonstrated at 
94 GHz at room temperature using n-type GaAs (94 GHz 
reflection beam isolator) . Then -tdLe final experiments 
were performed at the wave length of 337 pm using n- 
type InSb (IR reflection beam isolator) . 

8-1-1 GaAs 

O. Madelung [1964] quoted the experimental results 
for the effective mass of electrons in n-type GaAs as 
m*=0.067m^ and the static dielectric constant as 
Kj^=12.53. The mobility of the free carriers and the 
carrier concentration are' estimated by either Hall 
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measurement or van der Pauw measurement [van der Pauw, 
1958] . For nondegenerate n-type semiconductors the 
Hall coefficient is given by 


H ne 


( 8 - 1 ) 


where the factor r_ is determined by the scattering 

K 

mechanism. For n-type GaAs it is reasonable to assume 
that the polar optical scattering is dominant at room 
temperature. Since the relaxation time approximation 
fails for polar scattering, the calculation of r^^ for 
the case of polar scattering is not simple. B.F. Lewis 
and E.H. Sondheimer [1955] have calculated r^^ for this 
case using a variational method, and the factor r^^ for 
polar optical scattering is around 1.01 at room tempera- 
ture. Geometric effects due to the shape of the Hall 
sample and variation of Hall mobility with magnetic 
field are also discussed elsewhere [Madelung, 1964]'. 

The van der Pauw measurement, on the other hand, does 
not depend on the geometry of a planar sample or the 
position of the contacts . Hence the whole sample used 
in the experiment was used in the van der Pauw measure- 
ment. The parameters of n-type GaAs used in the experi- 
ment are summarized in Table [8-1] . 
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8-1-2 InSb 

The effective mass measured from the Faraday 
effects at the room temperature in n-type InSb was 
m*=0.021m^, and the value taken for static dielectric 
constant for n-type InSb is = 17.88 at room tempera- 
ture [Hogarth, 1965] . The mobility and carrier density 
of n-type InSb were measured by performing either Hall 
measurement or van der Pauw measurement [van der Pauw, 
1958] . It is known that in intrinsic InSb above 200°K 
both acoustic and polar optical mode scatterings are 
significant- The ratio of Hall mobility to the conduc- 
tion mobility r has been evaluated from a variational 
transport calculation by H. Ehrenreich [1959], and r^^ is 
found to be 1.01 at room temperature. After taking into 
account^this correction, the parameters of n-type InSb 
at the temperature of 284°K and 270. 6®K are listed in 
Table [8-1] . 

8-2 Experimental Set-up and Procedure 

The experimental apparatus is designed to measure 
the non-reciprocal reflection from the surface of a 
semiconductor CGaAs or InSb) covered by a dielectric 
medium. In the millimeter wavelength region, the 
experiments were performed at 94 GHz, and n-type GaAs 
was used as a plasma medium. In the submillimeter 
wavelength region, the experiments were performed at 



TABLE 8-1. Electric Properties of n-type GaAs and InSb 



GaAs 

InSb 

Temperature 

’300°K 

284°K 

270. 6“K 

Electron Effective 
Mass 

0.067m 

o 

0 . 021ra 

o 

0 . 021m 

o 

Dielectric 

Constant 

12.53 

17.88 

17.88 

Electron Mobility 

6 . lxl0^cm^/V~sec 

7 . 4xl0^cm^/V-sec 

7.9xl0^cm^/V-sec 

Carrier Density 

Id 

7.1xl0-^^/cm 

1.25xl0^®/cm^ 

1.06xl0^®/cm^ 
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337 um from HCN lasers, and n-type InSb was nsed as a 
plasma medium. 

8~2~1 94 GHz Reflection Beam Isolator 

The experimental apparatus shown in Fig. [8-1] is 
designed to measure the reflection from the interface 
between a dielectric and solid state plasma. The semi- 
conductor used for the experiment is n-type 1.44 £i-cm 
GaAs which was mechanically polished and chemically 
etched [H2SO^C4) + H^O^Cl) + H^Od)]. Pyramidi ally-cut 
dielectric (Stycast e=15e^, Emerson and Cuming) was 
used as a "lossless" dielectric medium, and the mirror- 
polished GaAs was mounted on the bottom of the pyramid 
as shown in Fig. [8-2] . Quarter-wave matching plates 
(Stycast, Emerson and Cuming) were placed on the incident 
and transmitted- sides . This device was placed at 300®K 
in the 1 inch air gap of a Varian 6- inch magnet which 
provides the magnetic field up to 15 kG. 

The pyramidal horns, whose beam angle is around 
7.6°, were used to provide the incident "plane" wave 
polarized in the plane of incidence, and as a detector 
pickup. These horns were set about 10 cm away from the 
sample in order to insure that the em wave incident on 
the surface of the dielectric was plane. a lock-in 
amplifier was used, to amplify and read the signal from 
the crystal detector. In order to reduce the effects of 




Fig. 8 t- 1. Block Diagram for 94 GHs EKperiroental Set-up. 











Fig. 8-2. 94 GHz Reflection Beam Isolator 
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stray signals, a window frame made out of Al foil was 
placed on the dielectric pyramid at the side of inci- 
dence . 

The reflected power was measured as a function of 
magnetic field at a fixed angle. The measurements were 
repeated in 5 degree increments from 30® to 75®. Total 
insertion loss due to the isolator was somewhat diffi- 
cult to measure because of use of pyramidal horns. The 
GaAs was replaced by a metallic conductor (polished 
brass) , and the reflected power was measured again. 

The net insertion loss from absorption due to the GaAs 
was then determined by comparing the power reflected 
from GaAs and that reflected from the metallic conductor 

8-2^2 IR Reflection Beam Isolator 

The block diagram of the experimental apparatus 
for the measurement of the far-infrared non-reciprocal 
reflection is shown in Fig. [8-3]. The InSb was 
mechanically polished and chemically etched. This semi- 
conductor was covered with a high density polyethylene 
layer with a thickness of 0.1 to 0.3 mm as shown in 
Fig. [8-4]. The sample was cooled by the thermoelectric 
cooler in order to increase the mobility and to decrease 
the carrier concentration of intrinsic InSb. One side 
of thermoelectric cooler was water-jacketed, and the 
lowest temperature obtained on its surface was -2.4®C. 




Fig. 8--3. Block Diagram for IR Experimental Set-up 
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Fig, 8”4. IR Reflection Beam Isolator 
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This devic^ was placed in the 1 inch airgap of a 6 inch 
Varian magnet which provided a maximum magnetic field 
of 15 kG. 

An HCN laser of 337 ym was constructed as an 
infrared source. The laser V7as operated primarily at 
the wavelength of 337 ym, with typical stable output 
of approximately 10 • mW CW. The linearly polarized out- 
put was extracted from a beam splitter. A lens was 
used as a window in order to focus the IP. signal onto 
the sample, and the beam angle was approximately 3° , 

The reflected wave from the sample was received by a 
pyroelectric detector which was calibrated by a thermo- 
pile. The detected signal was amplified by a lock-in 
amplifier and then recorded on an x-y recorder as a 
function of applied magnetic field- A window frame made 
out of the lossy material (Eccosorb, Emerson and Cuming) 
was used in order to reduce stray IR signals. To 
reduce the effects of stray reflection, the area around 
the sample exposed to radiation was also covered with 
the lossy material (Eccosorb, Emerson and Cuming) . 

The reflection and insertion loss from the inter- 
face were measured in a similar manner as before. The 
accuracy of the measurements was limited mainly by the 
stability of the HCN laser output power. Invar mirror 
spacers and water-cooled cathode provided the stable 
operation. The laser tube was cooled by air, and entire 
HCN laser system was placed on a massive alxominum 
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optical table. The variation of the output power from 
the HCN laser was under 5% over long periods. 

8-3 HCN Laser 

Far-infrared laser action produced by discharges 
through various gases containing carbon nitrogen- and 
hydrogen in their elemental forms have been extensively 
studied [Gebbie, et al . , 1964; Mathias, et al., 1965; 
Muller and Flesher, 1967] . This chapter presents 
briefly the design, construction and operation of a 
stable HCN laser, and the details are given elsewhere 
[Kanda, 1971] . 

8-3-1 Physical Mechanism of HCN Laser 

Several explanations for the physical mechanisms 
of HCN laser have been proposed on the basis that the 
stimulated emissions at the wavelengths of 337 ym and 
311 ym are attributed to the transitions within rotat- 
ional and vibrational states of the HCN molecule. 

D.R. Lide and A.G. Maki [1967] have found a set of 
vibrational levels in HCN which provides a satisfactory 
explanation for the stimulated emission. The intense 
lines at 337 ym and 311 ym are explained as transitions 
involving the 11^0 and 04^0 vibrational states, which 
are mixed by a Coriolis perturbation. 
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8-3-2 Design of Resonator 

'A wide variety of combinations of mirror radii can 
be used with any gas laser discharge tube to provide 
modes having desired values of diameter and angular 
spread- Practical resonator configurations which are 
most commonly used are confocal, sem.iconfocal and long 
radius configurations. Obviously not all combinations 
of mirror radii and separations are useable. Many of 
them originate diverging wave fronts that effectively 
focus all energy away from the mirror at the opposite 
end and result in very high diffraction loss. It is 
also known that lasers tend to operate in one or m.ore 
high order modes if the geometry permits them to do so. 
The only way to force a laser to operate in a uniphase 
wave front is to design a laser so that all modes 
except the lowest order mode have diffraction losses 
sufficiently high so that the operation is below 
threshold. 

The resonator employed for our HCN laser is the 
semiconfocal configuration which consists of one 
spherical m.irror and one flat mirror. The resultant 
mode has a relatively large diameter at the spherical 
mirror and a relatively small diameter at the flat 
mirror. The advantage of a semiconfocal configuration 
is relatively stable adjustment of the mirrors, and the 
disadvantage is the relatively inefficient pumping since 
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the cone-shaped mode intersects only one-third of the 
total volume of the laser resonator. 

The mode dimensions, which are defined as the radii 
for which the electric field falls to e ^ of its maximiom 
value at the surface of mirrors, for the semiconfocal 
resonator employed for our HCN laser are 

= 2.22 cm , 

^2 = 1.5 cm 

2 

When the Fresnel number N = vj- is chosen to be 1.6, 
average power loss for mode is 0.4% and that for 

TEMj^q mode is 7% . Thus the ratio of the losses of the 
two lowest modes is 17.5, and the aperture, radius a is 
3.28 cm . 


8- 3~ 3 Design of Beam Splitter 

Although there are several methods for extracting 
energy from far-infrared gas lasers., a hole-coupling 

t 

scheme is most commonly used because of the simplicity 
of its' construction. Usual hole-coupling schemes are, 
however, found to have the following disadvantages: 

1) Diffraction losses are quite high; 

2) Mode distortion (especially of the 
fundamental mode TEM^^) is large? 

There are ambiguities in designing the 
size of the coupling hole. 


3 ) 
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In order to eliminate these disadvantages, the energy 
extraction method employed is a dielectric beam split- 
ter with the Michelson interferometer configuration, 
which is a taut polyethylene m.embrane and is set at an 
angle of 45° to the laser axis. The thickness of the 
polyethylene is chosen in order to provide constructive 
interference of the beams reflected from each of its 
surfaces. This thickness is given by an odd multiple 
of 

t = I (n^ - 1/2) , (8-2) 

where A is the vacuum wavelength of the laser radiation 
and n is the refractive index of the polyethylene. For 
our HCN laser {A = 337 ym, n = 1.5) the thickness of the 
polyethylene membrane is 64 ym. 

Some of the advantages of this polyethylene beam 
splitter are: 

1) Diffraction losses are kept at a minimum 
because the full internal beam diameter 
is utilized; 

2) Mode distortion is minimized; 

3) The output beam is linearly polarized, 
which is especially a useful feature for 
our experiments ; 

4) The output coupling is easily varied by 

adjusting the side mirror shown in 

1 

Fig. [8-5]. 
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A disadvantage of this scheme is a reduction of Q of 
the resonator because of a very slight loss of the 
polyethylene . 

The power reflectivity of the polyethylene film 
at the wavelength of 337 ym is about 4% so that the 
maximum fraction of power coupled out of the HCN laser 
is about 16%. This can be increased even more, if the 
laser gain permits, by using a multiple parallel mem- 
brane with an appropriate spacing. 

8 - 3-4 Final Laser System 

Fig. [8-5] shows the schematic of an HCN laser 
system, operating 337 ym and 311 ym. The main optical 
resonator is inside the evacuated tube and is fully 
adjustable from the outside through the vacuvim seals.. 
The alignment of the mirrors was accomplished by using 
a He-Ne gas laser. The water-cooled cathodes were 
constructed from the copper plumbing items and were 
hollow in order to increase the discharge area and thus 
to provide a uniform discharge. Strong lasing action 
has been observed at 337 ym and 311 .ym with several 
chemical compounds [Gebbie, et al., 1964; Mathias, et 
al., 1965; Muller and Flesher, 1967]. A mixture of 
methane and ammonia was used for our HCN laser. The 
HCN laser requires a gas flow and exhaust system since 
the electric discharge dissociates the gas molecules. 
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and the resulting constituents are not reusable. The 
gas discharge operates at about 1.4 kV while the 
firing voltage is about 2.5 kV, but these voltages 
depend upon the pressure and. flow of the methane 
ammonia mixture. Final specifications for our HCN 
laser systemi are as follows: 


Length of laser tube 

2.1 m , 

Length of usable discharge 

1.7m, 

Inside diameter of laser tube 

7.5 cm , 

Curvature of mirror 

4 m , 

Pressure of methane ammonia 
mixture 

0 . 8~5 mmHg , 

Current limiting resistor 

1.8 kfl , 

Firing voltage 

2~3 kV , 

Operating voltage 

1.2-1. 6 kV , 

Operating current 

0.3-0. 8 A , 

Power output 

16 mW at 337 ym , 


5 mW at 311 ym . 


8-4 Infrared Radiation Detector 

The detectors used for infrared radiation may be 
divided into two distinct groups, thermal detectors and 
photodetectors. The thermal detectors use the heating 
effect of the radiation and depend on the number of the 
effective quanta of the radiation absorbed. The 
characteristics of the thermal detectors are that their 
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performance is independent of wavelength, and that the 
time constant is relatively long. The photodetectors, 
on the other hand, make use of the quantum photoelectric 
effects and depend on the total energy absorbed. One of 
the great advantages of the photodetectors over thermal 
detectors is that much shorter time constants are ob- 
tainable . 

The infrared radiation detector employed in our 
experiments is the pyroelectric detector [Beerman, 1967]. 
The actual sensor is a pyroelectric crystal exhibiting 
spontaneous polarization which is temperature dependent. 
Although such materials as triglycine sulfate (TGS) , 
triglycine fluoberyllate (TGFB) , lithium sulfate, 

Rochelle salt, barium titanate and lead zirconate 
titanates (PZT) are all pyroelectric, TGS has given 
better performance than these other crystals. 

The pyroelectric detector used was a piece of TGS 

2 

4 mm in area and O'. 5 mm thick. The tv/o opposite sides 
of the rectangular shaped sample were silver coated to 
create a capacitance with a thermal sensitive dielectric 
constant. Electrodes are attached to the silver coat- 
ing and connected to the high input impedance 2 stage 
PET amplifier. 

The pyroelectric detector is less sensitive than 
photodetectors and has a long time constant. However, 
it can detect long wavelength radiation without 
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cryogenic cooling. Also the pyroelectric detector, 
which is in principle a pure capacitor, does not show. 
Johnson noise and is generally . limited by the ampli- 
fier noise. 



Chapter 9 


THEOEETICAL AND EXPERIMENTAL RESULTS 

This chapter presents the theoretical and experi- 
mental results- Two geometries were considered; GaAs- 
dielectric interface at 94 GHz; free space-dielectric 
coated InSb interface at 337 ym. In all cases the 
polarization of the incident electric field is in the 
plane of incidence, and the direction of propagation 
is transverse to the magnetic field. A discussion of 
the theoretical and experimental results and sources of 
experimental error are also given. 

9-1 94 GHz Reflection Beam Isolator 

Theoretical reflection loss from the interface 
between free space and GaAs as shown in Fig. [9-1] 
indicates that there is very small isolation with high 
insertion loss, and even for other reasonable choices 
of parameters n, y, B and 0 [Seaman, 1969]. Therefore 
this device does not’ look promising for device applica- 
tion. It was thought that the high lattice dielectric 
constant (Kj^=12.53 for GaAs and Kj^=17.88 for InSb) 
might have washed out a large non-reciprocal behavior 
observed in the ionosphere by N.F. Barber and D.D. 
Cromhie [1959]. Thus, the reflection from the inter- 
face between a dielectric medium and GaAs was studied 
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Fig. 9-1. Theoretical Reflection Loss of GaAs 
at 94 GHz as a Function of Incident 
Angle [Parameters of GaAs are listed 
in Table 8-1. Geometry of 94 GHz 
Reflection Isolator is shown in 
Fig. 7-la with B=13.2 kG] . 
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to achieve better isolation with smaller insertion 
loss . 

Theoretical and experimental reflection loss from 
GaAs with a semi-infinite dielectric ‘(K^=15) for the 
electric field polarized in the plane of incidence at 
94 GHz and at 300 °K are shown in Fig. [9-2] for a dc 
magnetic field of 13.2 kG as a function of incident 
angle, and in Fig. [9-3] at a fixed incident angle of 
57.5° as a function of a dc magnetic field. 

Isolation of 12 db with 11 db insertion loss was 
observed for the signal at an incident angle of 57.5° 
whereas the theory predicts 44 db isolation with 5 db 
insertion loss at an incident angle of 65°. The devia- 
tion between theoretical and experimental results were 
great. In the next several paragraphs/ possible sources 
of experimental error are discussed. 

The surface of GaAs was mechanically polished and 
was flat and smooth to within a small fraction of the 
wavelength at 94 GHz. The thickness of GaAs was 4.4 mm 
which was more than 9 skin depths , where the skin depth 
is defined as the depth at which the magnitude of a 
field quantity decreases to e ^ (about 36.9%) of its 
value at the surface. Therefore the reflection from 
the back side of the GaAs was negligible. 

The parameters of GaAs were determined by the 
van der Pauw measurement which should give an overall 
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Fig. 9-2. Theoretical and Experimental Reflection Loss 
of GaAs at 94 GHz as a Function of Incident 
Angle [Parameters of GaAs are listed in . 
Table 8-1. Geometry of 94 GHz Reflection 
■ Isolator is shown in Fig. 8-2 with B=13.2kG] 
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Fig. 9-3. Theoretical and Experimental Reflec 
94 GHz as a Function of Magnetic Fi 
GaAs are listed in Table 8-1. . Geom 
Reflection Isolator is shown in Fig 
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average value of the mobility and the carrier concen- 
tration of free carriers for the whole sample of GaAs . 
However, the measurements made on the various pieces 
indicated that the GaAs used in our experiment was not 
uniform and variations of 10% in these parameters were 
found to exist. This caused some troiible in the experi- 
ments . 

The 94 GHz single from the pyramidal horns used 
for transmitting and receiving has a substantial angular 
spread between 2® to 5° depending on the tilt of the 
sample to the horn, and therefore the reflection 
obtained from the experiments was an "average" value of 
reflection for an angular spread. The spread in inci- 
dent angles not only smeared out the sharpness of the 
reflection curve, but also reduced its peak as predicted 
from the theory (see Fig. [9-2]). 

The effects of stray reflections and transmissions 
were minimized by the Al window frame placed on the 
pyramidally- cut dielectric at the side of incidence. 

The size of the window was, however, only about 5>^7 
wavelengths , and the non-plane -wave behavior due to 
diffraction caused by this "small" window have washed 
out the sharpness of the reflection loss. This was 
also confirmed from the fact that a much smaller 
window, 3 wavelength square, completely destroyed the 
non- reciprocal phenomena . 
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In order to estimate the net insertion loss from 
absorption due to GaAs / the difference of reflected 
power from GaAs and from a metallic conductor was 
measured. Although the sample holder was made in such 
a way that it was possible to replace the GaAs by a 
metallic conductor without moving the sample holder 
itself, it was still difficult to position a metallic 
conductor exactly in the same position as was GaAs. 

This caused some error in estimating the incident power 
on the GaAs and thus the insertion loss. 

The largest reflection loss was observed at the 
incident angle of 57.5° whereas the theory predicted a 
peak at 65°. Since the pyramidal horn has a beam angle 
of 7.6°, it is quite difficult to determine a definite 
angle of- incidence, and thus the measurement of the 
incident angle of the wave from the normal was not 
accurate. 

9-2 IR Reflection Beam Isolator 

The primary experiment was conducted on the reflec- 
tion from the interface between free space and the InSb, 
and then the further experiments were carried on the 
reflection from the interface betv;een free space and a 
dielectric coated InSb. All theoretical calculations 
and experiments were performed at a wavelength of 337 um 
from HCN laser. The theoretical and experimental 



reflection loss from the interface between free space 
and InSb at 284°K are shown in Fig. [9-4] for a dc 
magnetic field of 15 kG as a function of incident angle, 
and in Fig. [9-5] at a fixed incident angle of 64® as 
a function of the dc magnetic field. General agreement 
between theoretical and experimental results is quite 
good although, as expected, these results indicate 
that this device does not look promising for the device 
application. It was thought that the high lattice 
constant in InSb (Kj^=17.88) have washed out the large 
non-reciprocal reflection, and therefore the reflection 
from the interface between a semi-infinite dielectric 
medium and a solid state plasma was considered just as 
for the 94 GHz reflection beam isolator. The theoreti- 
cal reflection loss from the interface between a 
dielectric layer (Kj^=30) and InSb at a wavelength of 
337 ym and at 300®K with a dc magnetic field of 13.2 kG 
was calculated using the theory given in Section 7-2 
and is shown as a function of incident angle in Fig. 
[9-6]. The isolation of 34 db with 0.5 db insertion 
loss is predicted at an incident angle of 85.5°. The 
necessity of the dielectric medium whose dielectric 
constant is 30, however, makes this experiment difficult. 

Rutile, [Ti02], barium titanate [BaTiO^] , barium- 

1 

strontium titanate [ (Ba-Sr) Tio^] and a large number of 
bismuth compounds are known to have a high dielectric 
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Fig. 9-4. Theoretical and Experimental Reflection Loss of InSb at 
337 ym as a Function of Incident Angle [Parameters of 
InSb (T=284°K) are listed in Table 8-1. Geometry of IR 
Reflection Isolator is shown in Fig. 7-la, with_R=15kG] 
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Theoretical and Experimental Reflection Loss of InSb at 
337 ym as a Function of Magnetic Field [Parameters of 
InSb (T=284°K) are listed in Table 8-1. Geometry of IR 
Reflection Isolator is shown in Fig. 7-la with 0=64°]. 
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Fig. 9-6. Theoretical Reflection Loss of InSb at 
337 um as a Function of Incident Angles 
[Parameters of InSb (T=284°K] are 
listed' in Table 8-1. Geometry of IR 
Reflection Isolator is shown in 
Fig. 7-lb with Kj^=30 and B=13.2kG]. 
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constant. But none of these high dielectric materials 
are transparent at the wavelength of 337 ym. The few 
materials, known to be transparent at this wavelength, 
(some semiconductors cind some dielectric materials such 
as polyethylene, teflon and quartz) have a typical 
dielectric constant of around 10 or less, and none of 
these is even close to what is required for this experi- 
ment (Kj^=30) . 

In order to overcome this difficulty, a thin 
dielectric layer placed on top of InSb was introduced 
as a matching transformer. Thus the reflection from 
the interface between free space and dielectric coated 
InSb was studied. The idea of this configuration was 
to create the same effect on the reflection from the 
interface between a high dielectric medium (K^=30) and 
InSb by adjusting the thickness of a dielectric layer, 
whose dielectric constant is around 2, placed on top of 
InSb. The dielectric material chosen for this layer 
was high density polyethylene whose dielectric constant 
is 2.27, and its loss tangent (tan S) is as low as 
1.3x10 ^ at the wavelength of 337 ym [Breeden, 1968]. 

The possibilities of using a thin layer with a different 
dielectric constant for the purpose of matching have 
been considered and are given in Section 10-2. 

The theoretical and experimental reflection loss 
from the interface betv^een free space and dielectric 
coated InSb at 284°K are shown in Fig. [9-7] for the dc 
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magnetic field of 15 kG as a function of incident 
angle, and in Fig. [9-8] at a fixed angle of 65° as a 
function of dc magnetic field. The thickness of the 
high density polyethylene assimed in the theoretical 
calculations and used in the experiments was 250 pm. 

The general shapes of the curves for the theoretical 
and experimental results agree fairly well except at 
large incident angles. Since at large incident angles 
only a small IR signal was incident on InSb and con- 
sequently only a small IR signal was reflected off from 
InSb, the experimental data at incident angles over 75° 
were questionable. Therefore effort was made to shift 
the isolation peak to a lower angle, as discussed in the 
next paragraph. 

F rom the heuristic argument the isolation peak 
takes place when the impedance of the dielectric medium 
is matched with that of the InSb. The impedance of the 
dielectric medium is given by geos 6 where g is the 
intrinsic impedance of the 'dielectric medium and 0 is the 
incident angle. Similarly the impedance of the InSb is 
given by g^ cos0^ where g^ is its intrinsic impedance 
and 0^ is the angle of refraction in the InSb. When 
the intrinsic impedance g of the dielectric medium and 
the angle of refraction 0^ in the InSb are given, the 
incident angle 6 can be shifted to a lower angle by 
increasing the intrinsic impedance g^ of the InSb. 
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Theoretical and Experimental Reflec- 
tion Loss of InSb at 337 ym as a 
Function of Incident Angle [Parameters 
of InSb (T=284°K) are listed in Table 
8-1. Geometry of IP. Reflection 
Isolator is shown in Fig. 8-4 with 
A=250ym and B=15kG] . 
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Fig. 9-8. Theoretical and Experimental Reflection 
Loss for InSb at 337 ym as a Function 
of Magnetic Field [Parameters of InSb 
(T=284°K) are listed in Table 8-1. 
Geometry of IR Reflection Isolator is 
shown in Pig. 8-4 with ii=250ym and 
e=65°] . 
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It is found theoretically that the intrinsic impedance 
of the InSb with a dc magnetic field can be increased 
with an increase of its mobility and a decrease of its 
carrier concentration. It is, therefore, expected that 
the peak of the isolation curve can be shifted to a 
lower incident angle by few degrees by increasing the 
mobility and decreasing the carrier concentration. 

One way to achieve this is to cool the InSb with a 
thermoelectric cooler. By cooling the InSb from 284°K 
to 270. 6°K, the mobility increases by 7% and the carrier 
concentration decreases by 15%. This, in turn, shifts 
the peak of the isolation curve by 2®, which turns out 
to be a too small amount to ease experimental difficul- 
ties. In fact, a drastic shift of the isolation curve 
to a lower incident angle cannot be achieved only by 
increasing the mobility and decreasing the carrier 
concentration so long as polyethylene {¥.^ 2 , 21 ) is used 
as a matching transformer; further considerations are 
made in Section 10-2. 

The theoretical and experimental reflection loss 
from the interface between free space and a dielectric 
coated InSb at 270. 6®K is shown in Fig. [9-9] for a dc 
magnetic field of 12 kG as a function of incident 
angle, and in Fig. [9-10]' for a fixed incident angle of 
60° as a function of dc magnetic field. The thickness 
of the high density polyethylene assigned in the 
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Fig* 9-9. Theoretical and Experimental Reflection 
Loss for InSb at 337 um as a Function of 
Incident Angle [Parameters of InSb 
(T=270.6°K) are listed in Table 8-1. 
Geometry of IR Reflection Isolator is 
shown in Fig. 8-4 with Jl= 125 yin' and 
B=12kG] . 
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Fig. 9-10. Theoretical and Experimental Reflection Loss for InSb 
at 337 ym as a Function of Magnetic Field- [Parameters 
of InSb (T=270.6°K) are listed in Table 8-1. 

Geometry of IR Reflection Isolator is shov/n in 
Fig. 8-4 with £=125ym and 0=60'’]. 
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theoretical calculations and in the experiments was 
125 ym. The best data obtained experimentally as an 
application to an isolator were 4.2 db isolation with 
2 db insertion loss at an incident angle of 77.5® with 
the dc magnetic field of 15 kG. With this magnetic 
field, the theory predicted 16 db isolation with 1 db 
insertion loss at an incident angle of 86°. Large 
discrepancies between theoretical and experimental 
results take place at incident angles over 80°, as 
discussed earlier. Therefore, further considerations 
are made in Section 10-2 in order to shift the peak cf 
the isolation curve to a lower incident angle. Except 
at incident angles over 80°, the experiments seem to 
indicate that the theory is valid, and general agree- 
ment between theoretical and experimental results is 
quite good within experimental errors . In the next 
several paragraphs, possible sources of experimental 
error are discussed. 

InSb was flat and smooth to within a small fraction 
of the wavelength (337 ym) . Since the thickness of the 
InSb was 1 . 4 mm which was approximately equivalent to 
70 skin depths, reflection from the back side of the 
InSb was negligible. The parameters of the InSb were 
measured by the van der Pauv; method which gives accurate 
average values for whole samples. The laser output at 
337 ]xm was carefully separated from that at 311 ym by 
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adjusting the fine tuning of the resonator. The beam 
splitter was used to extract the polarized output 
power, and the degree of polarization was ascertained 
to be over 99.9% [Evenson, 1969] . The long time stab- 
ility of the output power of HCN laser was reduced to 
under 5%. The errors introduced in the results from 
these sources are, therefore, quite negligible. 

The output power from HCN laser has an angular 
spread of around 1° depending on the tilt of the sample 
to the window of the laser. Therefore the reflection 
obtained from the experiments V7as an "average" value of 
reflection for an angular spread. The spread in 
incident angles not only smeared out the sharpness, but 
also reduced the peak of the reflection curve. There- 
fore the- error caused by this source was significant 
especially at large incident angles where the reflection 
coefficient changes rapidly with incident angles (see 
Fig. [9-7] and [9-9]). 

The thin layer of high density polyethylene placed 
on top of the InSb was used as an impedance transformer. 
Since the wavelength was quite short (337 ym) and 
interesting experiments had to be performed at large 
incident angles, the accuracy of the thickness of the 
high density polyethylene is important. Since the 
polyethylene was cut by milling, some error was expected 
due to the limit of machining accuracy. 
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In order to eliminate the stray reflections and 
transmission, a window frame made out of lossy material 
(Eccosorb, Emerson and Cuming) was placed in front of 
the InSb. This ensured that the detected signal was 
indeed the signal reflected from the InSb, and not 
stray transmission or reflection aside from InSb, 

The semiconductor was replaced by a metallic 
conductor (polished brass) , and the reflected power 
was measured again. The net insertion loss from 
absorption due to the InSb was then determined by 
comparing the reflected IR signal from the InSb and 
that from the metal. If the conductor was not positioned 
exactly in the same place as was InSb, the IR signal 
incident on each surface would be different, thus causing 
some error in estimating the insertion loss . 



chapter 10 


DISCUSSIONS AbTD CONCLUSIONS OF PART II 

In this chapter the behavior of the em waves 
reflected from a magnetoplasma is discussed. The 
physical quantities, such as the ellipticity of the 
electric field and the orbit of electrons in a magneto- 
plasma excited at the plasma surface due to an incident 
em waves are evaluated. In all cases the polarization 
of the incident electric field is in the plane of inci- 
dence, and the direction of propagation is transverse 
to a dc magnetic field. It is hoped that the evaluation 
of these physical quantities will lead to the physical 
understanding of the non-reciprocal phenomena. Finally 
applications of the device and possible refinements for 
isolators are also discussed. 

10-1 Physical Explanation of Non-reciprocity 

A dc magnetic field makes the plasma anisotropic. 
The origin of this anisotropy is due to the motion of 
the electrons in a dc magnetic field. A component of 
electric field perpendicular to a dc magnetic field can 
give rise to the transverse motion of electrons. This 
transverse velocity interacts with a dc magnetic field 
to produce motion in the Hall or vxB direction. As a 
result, off-diagonal terms appear in the matrix relating 
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the velocity (and the current) to an electric field. 
Therefore conductivity and dielectric constant become 
tensors . It can be shown that the nature of the 
magnetic force is such that the dielectric tensor is 
antisymmetric, and off-diagonal terms contain terms 
linear in the dc magnetic field as long as the field is 
small. This anisotropic effect in the plasma is one of 
the essential factors required for non-reciprocity. 

As discussed in Section 2-3, in the case of propaga- 
tion transverse to a dc magnetic field, the propagation 
constant is reciprocal. That is, reversal of the direc- - 
tion of propagation, or alternatively the direction of 
a dc magnetic field, does not change the propagation 
constant. Thus the question arises why non-reciprocal 
reflect ion takes place. In order to settle this, let 
us consider, first of all, the excitation of em v/aves 
due to an incident plane wave on a magnetoplasma. 

The non-reciprocity arises from the interface 
between free space and the plasma. The coupling between 
the X and y components of the electric field through 
the dc magnetic field is expected from the Hall current 
due to the anisotropic property of the tensor dielectric 
constant. Since the Plall current has a definite sense 
with respect to the dc magnetic field, the amplitude 
of the electric field contributed by the Hall current 
contains the term linear in the field. That is, for 
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one direction of magnetic field the x component of 
electric field, for example, is coupled to the y com- 
ponent of electric field through a positive coefficient, 
and for the other direction of magnetic field is coupled 
through a negative coefficient as shown in Fig. [10-1], 
Therefore, total electric field in the plasma is not 
reciprocal when the dc magnetic field is reversed. 
Consequently, the reflection coefficient at the inter- 
face between free space and InSb is non-reciprocal. 

An explanation for the non- reciprocal reflection 
based on the changes in electron orbits in a magneto- 
plasma was proposed by K. Davies [1969] , His argument 
is extended here in order to give a more convincing 
physical explanation for the non- reciprocal reflection. 
Consider' the situation v/hen two identical waves are 
traveling through a plasma medi\im at equal but opposite 
angles to the vertical (see Fig. [7-la]) . In the absence 
of a dc magnetic field, the ratio of the x and y 
components of the electric field in a plasma is given 
by 

• ^ ^ -_sine 

^x . 2 , 

(K I - sin 0 ) 

where 

2 
CO 

K, = V 

]_ L w(w-3V) 


( 10 - 2 ) 
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Fig. 10-1. Electric Field Excited in 
a Magnetoplasma (Electric 
fields excited through the 
Hall current are indicated 
by *■ and — — ) . 
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Thus if the plasma is lossy, i.e. is complex, the 
electric field in the plasma is elliptically polari'zed 
in the plane of incidence as shown in Fig. [10-2] , 

(note the difference in x and y scales) . Hence the 
orbits of electrons responding to this electric field 
are also elliptic. These elliptical motions can be 
decomposed into two oppositely rotating circular orbits 
with different radii. When a transverse dc magnetic 
field is applied, the radius of each circle is increased 
or decreased depending on the direction of a transverse 
dc magnetic field. The resulting electron orbits ob- 
tained by summing the two oppositely rotating circular 
orbits may be drastically different from those without 
3 transverse dc magnetic field. That is, for one sense 
of magnetic field, the ellipse will be flatter than the 
case with no magnetic field, and for the other direction 
of magnetic field the ellipse will be more round, as 
indicated pictorially in Fig. [10-3a] . 

Fig. [10-4] shows the validity of the above 
physical interpretation. Fig. [10-4a] shows the orbits 
of electrons in a magnetoplasma with a dc magnetic field 
of 2 kG when non-reciprocity is quite small. It is ' ' 
found in Fig. [10-4a] that the orbits of electrons in 
a magnetoplasma are almost reciprocal with reversal of 
the dc magnetic field and hence small non-reciprocal 
reflection is expected. Fig. [10-4b] shows the orbits 
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Pig. 10-3. Pictorial Explanation for 
Non-reciprocal Reflection. 




Fig. 10-4. Orbits of Electrons in the InSb W3 th 

Reversal of dc Magnetic Field fParameters 
of InSb (T=284°K) are listed in Table 8-1 
Geometry considered is shown in Fig. 7-la 
with 0=50°]. 
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of electrons in a m^gnetoplasma with the dc magnetic 
field qf 15 kG when no'n- reciprocity is quite large. 

Fig. [10-4b] indicates large change in electron orbits 
with reversal of the dc magnetic field, and therefore 
large non-reciprocal reflection is observed. 

The explanation based on changes in electron orbits 
with reversal of the magnetic field also explains why 
collisions must be present,* if not, the reversal of the 
direction of propagation does not change the magnitude 
of, the reflection coefficient but only its phase-. Also 
if the collision rate is too high, only small non- 
reciprocal reflection takes place. The orbit of elect- 
rons in a collision-less medium or in a semiconductor 
with a small mobility is almost linear, as already 
shown in Fig. [10-2] , (note difference in x and y 
scales) . This linear polarization can be decomposed 
in right and left circular polarizations \vith equal 
radii. Application of a dc magnetic field causes one 
polarization to increase in magnitude and the other 
to decrease, and hence the polarization becomes 
elliptic. Reversing the field causes the whole process 
to reverse, but elliptic orbits of an electron for two 
directions of the dc magnetic field are identical. 
Therefore no non-reciprocal refledtion is expected, but 
there is a change in phase, as already shown in Fig. 


• [10-3b]. 
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From the discussions given above, it appears that 
a maximum difference in a reflection from a magneto- 
plasma with reversal of the dc magnetic field would take 
place if the orbit of electrons with B=0 is nearly • 
circular. Of course, it cannot be completely circular 
because of loss in the plasma as shown in Fig. [10-2] . 
Fig. [10-5a] shov/s that the orbit of electrons in the 
InSb with a lov? motility of 50 cm^/V— sec, where the 
reflection is found to be nearly reciprocal, is almost 
linear. On the other hand. Pig, [10— 5b] shows that 
the orbit of electrons in the InSb v/ith a much higher 
mobility of 7 10 ’cm^/V-sec , where an isolation of 
1 db is observed with B=15kG at an incident ancle of 
60°, is elliptic. The isolation given above is quite 
small because only the simplest case, reflection from 
the interface between free space and InSb, was 
considered. Large isolation v?as , however, observed 
in the reflection from the interface between free 
space and dielectric coated InSb as given in Section 
9-2. 


Thus, let us investigate the eccentricity e of 
the orbit of electrons v/ith B=0 . (Recall that e=l 
for a straight line, e<l for an ellipse and e=0 for 
a circle) , 


. For our case the eccentricity e is given by 

2 2 'll ^ 


1 - 


1 - 1 ^. 


1 - M 


_ iKlcos^ ^ (l+„2) /2M:COSl)2^ 

- M V X - 


1 - 


(10-3) 
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Fig. 10-5. Change in Electron Orbits in InSb 

with no Magnetic Field for Different 
Mobilities [Other parameters of InSb 
(T=284°) are listed in Table 8-1. 
Geometry considered is shown in 
Fig. 7-la with 0=50° and B=0 1 . 
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where 
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The eccentricity e as a function of incident angle 9 
is shown in Fig. [10-6] for several boundary conditions. 
Since the eccentricity decreases, i.e. the orbit of 
electrons in a plasma becomes more circular, with 
increasin'^ incident angle, it is expected that the 
difference in reflection from a magnetoplasma with 
reversal of the dc magnetic field should increase with 
incident angle. Kov/ever, at large incident angles 
(over 85°) the reflection is quite large and hence only 
a small em wave is transm.itted in a plasma. Thus the 
change in reflection v/ith reversal of the dc magnetic 
field decreases at large incident angles simply 
because the wave does not interact much with the 
plasma. It is also found that the eccentricities of 
the orbits of electrons for the cases of dielectric- 
plasma and free space-dielectric coated plasma 
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INCIDENT ANGLE (degrees) 


Fig. 10-6. Eccentricity e of Electron Orbits 
in InSb as a Function of Incident 
Angle [Parameters of InSb (T=284°K) 
are listed in Table 8-1. Geometry 
considered is shown in Fig. 7-1 
with 0=50° and B=0] . 



125 


interfaces are much smaller than that for the case of 
free space-plasma interface as shown in Fig. [10-6] . 

That is, the orbit of the electrons is much less 
linear, and hence greatly enhanced non-reciprocal 
reflection can be observed for the cases of dielectric- 
plasma and free space-dielectric coated plasma inter- 
faces . Therefore, the combined effects of both 
collisions, which are required to have an elliptic 
motion of electrons in a plasma without a magnetic 
field, and a dc magnetic field produce non-reciprocal 
effects on the electron orbits when the dc magnetic 
field is reversed, and thus explain the non-reciprocal 
reflection. 

It should be pointed out that, since the electric 
field i's~elliptically polarized in a lossy magneto- 
olasma, there is a component of the electric field along 
the direction of the propagation- Because of this 
longitudinal component, the existence of space charge 
might be expected. The deviation in carrier concentra- 
tion from equilibrium can be calculated from the 
Poisson’s equation and is given by 

, wy £ 0 ^ 

n = K — ■ °-S . sin0 Im [ (^ --sin 0)] (10-6) 

o e A 

This turns out to be around 10 /cm whereas the equili- 

16 3 

brium carrier concentration is 10 /cm . Therefore, the 
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charge neutrality conc3ition holds. This is also confirm- 
ed from the fact that tlie dielectric relaxation time is 
on the order of 10 to 10 sec, on the other hand, 

the period of the phenomena that we are interested in is 

-1 1 ’ ' - . 
around 10 " to 10 ' ' sec. Thus the space charge can be 

neglected, 

t 

10-2 Possible Refinements for Isolators 

Large discrepancies between theoretical and experi- 
mental restilts take place at incident angles over 80° 
due to experimental difficulties as discussed in Section 
9-2. Thus further considerations are made in this 
section in order to shift the peak of the isolation 
curve to a lower incident angle. 

First making use of a thin dielectric layer with a 
dielectric constant less than unity or even with a 
negative dielectric constant was considered in order to 
have no reflection at a reasonable incident angle. 

Such a layer can be realized by using either a semi- 
conductor or an artificial dielectric. For a semi- 
conductor, the carrier concentration must be chosen so 
that total relative dielectric constant contributed 
from lattice as well as free carriers is less than unity 
or even negative. For the ideal case when t=“ and 
Eq. (2-14) , (2-15) and (2-16) become 
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(t) 




p 


< 1 


w 


(10-7) 


= 0 , ( 10 - 8 ) 

if the carrier concentration of the semiconductor is 
chosen properly. A dielectric medium whose dielectric 
constant is less than unity can be also realized using 
an artificial dielectric. The arrays of conducting 
rods or plates were proposed as an artificial dielec- 
tric [Harvey, 1970]. When the dielectric layer on top 
of the InSb has a dielectric constant of 0.697, it is 
possible to create the condition that no reflection 
takes place at a reasonable incident angle (0=60®) as 
shown in Fig. [10-7], However, the peak of the isola- 
tion curve— is very critical with incident angle, and 
hence with a carrier density and mobility, and with dc 
magnetic field. This sharpness of the isolation curve 
is essentially due to the critical impedance match 
which results from the large impedance difference 
between free space and the InSb. However, the sharp- 
ness of the isolation curve could be made broader by 
the use of multiple layer matching transformers. By 
tolerating a certain mismatch at the central incident 
angle, it is possible to match at N incident angles by 
using N matching transformers to give an even wider 
bandwidth [Young, 1959]. It may be possible to make 



REFLECTION LOSS (db) 



Fig. 10-7. Theoretical Reflection Loss of 
InSb at 337 Tim as a Function of 
Incident Angle [Parameters of 
InSb (T=289°) are listed in 
Table 8-1. Geometry of IR 
Reflection Isolator is shown 
in Fig. 8-4 vrith B=16kG, 

K =0.697 and Jl=2 76ym] . 
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use of this sharpness of the isolation for the develop- 
ment of the angle resolving devices. 

The utilization of Alfven effect was also consider- 
ed. The typical ' materials considered here are bismuth 
and bismuth- antimony alloy which have a relatively high 
mobility, and low carrier concentration for both elect- 
rons and holes . The existence of both electrons and 
holes with almost equal mobilities and masses requires 
some modifications in the calculation of the reflection 
coefficient. This might, in turn, help in shifting the 
peak of the isolation curve to a lower incident angle - 
It turns out, unfortunately, that the cancellation takes 
place in the off-diagonal component of the relative 
dielectric tensor, K^, which plays a vital role in the 
non- reciprocal phenomena. Therefore the relative di- 
electric tensor becomes diagonal, and non-reciprocal 
phenomena cannot be expected. 

10-3 Summary 

The theoretical and experimental investigations 
reported in Part II were performed for the purpose of 
the development of the reflection beam isolator at 
millimeter and submillimeter wavelengths using non- 
reciprocal reflection of em waves incident on a semi- 
conductor. The best data obtained in our experiments 
at 94 GHz were 12 db isolation with 11 db insertion 
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loss, and corresponding theoretical data were 44 db 
isolation with 5 db insertion loss. At a wavelength 
of 337 yin the experiments demonstrated 4.1 db isolation 
with 2 db insertion loss, whereas the theory predicts 
16 db isolation with an insertion loss of less than 
1 db. However, since theoretical values given above 
occur at very large incident angles, difficulties arose 
in the experim.ents - Otherwise, the experimental results 
indicated the validity of the theory in general. The 
theoretical result for a wavelength of 118 ym shows a 
very broad reflection loss curve (about 4 db) as a 
function of incident angle and an excellent ratio of 
isolation to insertion loss of about 50 in db scale as 
shovm in Fig. [10-8]. The reverse loss can be increased 
by the multiple reflection. This result indicates that 
the development of the practical reflection beam isola- 
tor should be possible using solid state magnetoplasmas. 

One of the great advantages of utilization of semi- 
conductors in the developm.ent of non-reciprocal devices 
lies in the flexibility of design that comes with a 
wide range of semiconductor parameters. Thus it is 
possible to expect significant improvement of the non- 
reciprocal semiconductor devices as a result of continued 
research and development effort. 
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INCIDENT ANGLE (degrees) 

Fig. 10-8. Theoretical Reflection Loss for InSb at llSym as a Function 
of Incident Angle [Parameter of InSb (T=284°K) are listed 
in Table 8-1. Geometry of IR Reflection Isolator is shown 
in Fig. 7-lb with Kj^=50 and B=13.2kG]. 


131 



APPENDIX A. ENERGY DEPENDENCE OP RELAXATION TIME • 


It is known that the relaxation time is dependent 
on the energy of ' electrons . This Appendix will esti- 
mate the error caused by the assumption that the 
relaxation time is not dependent on the energy of 
electrons . 

The equations given in this thesis, such as propa- 
gation constants, reflection coefficients and etc., can 
be expressed in terms of the elements of the relative 

When the relaxation 
time is a function of energy, the Boltzmann equation 
must be used to derive the relative dielectric tensor. 
The solution of the Boltzmann equation requires the 
averaging of the quantities containing the relaxation 
time in the relative dielectric tensor with respect to 
energy. But otherwise, the forms of the elements of 
the relative dielectric tensor are the same as for the 
case when the relaxation time is not a function of , 
energy. The results are 


aieiectrxc tensor, K 


1 ' 


and K 


K, = 


1 


K - i-P 

L to 


( jtOT + l) T 


( jtOT+1) ^ + (w^t) ^ 


> , (A-1) 


= 




to T 
c 


(jtOT + l)'^ + 


(A-2) 
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and 


K, 


= - 




(0 


"j COT+1' 


(A- 3) 


where <> indicates the appropriate averaging of x with 
respect to energy. Since the averaging of x in these 
expressions is quite complicated, we consider these in 
the following cases : 


1. High Frequency Approximation (wx>>l) 

1.1 High Field Case (to >>to) : 

c 


Ki = 




to 

[1 + 4- <i>] 


to 


:)to X 


(A-4) 


K. 


. 2 

to to 


(A- 5) 


K I i = K 


. !!p_ [1 „ i_ <i>] 

L ^2 jco T ' 


(A- 6) 


1.2 Lov 7 Field Case (to>>to ) : 

c 


to 


Kj = [1-4— <^>] 

1 ^ ,.2 .X ' 


to 


(A- 7) 


K„ = 


. 2 

ItO to T 

[1 - 2_ ^1 
0)3 3^) X J 


(A- 8) 


K| 


to - 

= K_ ^ [1 - ™ <i>] 

L ^2 30 ) X J 


(A- 9) 
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2. Low Frequency Approximation (wt<< 1) 
2.1 High Field Case (w^>>w) : 




W T 
C 


(A-10) 


K.. = 


. 2 
-J(0 

P 


.X 0)W 


(A-11) 


. 2 

3^ 


K 1 I = K. - 


L W 


<T(1 - jWT)> 


(A-12) 


2.2 Low Field Case (co>>w ): 

c 


K| = Kj ^ 

11“ 


<T (1 - jUT) > 


(A-13) 


3 - 2 

a - 2jwT)> 


(A~14) 


K I I ■ = K 


3“ 


P 


L 0) 


<T (1 - jOJT)> 


(A-15) 


Since in our experiments , the value pf lOT was 
greater than unity {wT'“4) , for the first order approxi- 
mation, the values of Kj^, and Kj j for the case of 
0JT>>1 can be used in the development of the theories. 
If the relaxation time x is isotropic, then x may be 

expressed as a function of energy e , rather than as a 

* 

functxon of the wave vector k. For the non-degenerate 
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case/ the average value of the reciprocal of the relaxa- 
tion T, which occurs in the formula of K's are then 
evaluated by 

(A-16) 

0 


< 1 > 


1 <, 3/2 
— e e 

T 


e 

kT 


ds 


^ 3/2 
£ e 


£ 

kT 


de 


When the energy dependence of t is assumed to be 
-S 

T=ae , then 


<T> 


<i> 

T 



- s)r(| + s) 

5 ^ 

[r(|)] 


(A-17) 


For the acoustic scattering where S = j. 


<T> <— > 
T 


s= 


1 9it 


= 1.13 


(A-18) 


and for the ionized impurity scattering where S - - jr 


"""VI 3 = 17=3.38 . (A-19) 

2 

In the case of polar scattering, the average of <t><^> 
has not been calculated. Judging from the temperature 
dependence of the mobility of free carriers for n-type 
InSb and n-type GaAs at room temperature, the acoustic 
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scattering is probably dominant. Thus the maximum 
error caused by using <x> should.be 13% from Eq. (A-18) . 



APPENDIX B. CHARACTERISTIC EQUATION FOR CIRCULAR 

WAVEGUIDE CONTAINING ANNULAR PLASMA COLUMN 

Appendix B outlines the boundary value problem 

that results frora tlie configuration shown in Fia.(3-1) . 

In the region r<a, the wave equations describing the 

longitudinal components of em waves are 

vJe° + q^E° = 0 (B-1) 

+ q^H° = 0 (B-2) 

The appropriate solutions for the longitudinal field 
components for r<a are^ 

E° = UJ (qr) e^™*^ (B-3) 

z m 

and 

= VJ (qr)e^™‘^“^^ (B-4) 

z m ^ 

where U and V are arbitrary, complex constants. Correspond- 
ing transverse components of em waves are calculated from 
the equations : 


■±o 

^T 

(-r/q")V^E^ 4- 


(Bt5) 

- 

Ht - 

{-r/q2)V^H^ - 

(Jo3£^/q2) (J^xV^E^) 

(B-6) 


Then the transverse E field components for r<a are 

E° = [agEJ^(qr) - V j (qr) ] , (B-7) 

E° = [n|^EJj^{qr) + VCZoJ,;(qr)] . (b-8) 

Also the transverse H field components for r<a are 
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= * VqSJjljlqr)] (B-9) 


+. vj2jj^(qr)] 


jm<J)-rz ■ (B-10) 


In the region of the annular plasma column a<r<b. 
Maxwell ' s equations yield the coupled v/ave equations ; 


V^E + aE = bH„ 
T z z H 


(B-11) 


V„H + cH = dE 
T z z z 


(B-12) 


This coupled set of second-order equations can be trans- 


formed into uncoupled set of fourth-order equations: 


+ (a+c)vj + (ac-bd)] E^ = 0 
+ (a+c)v2 + (ac-bd)]H^ = 0 


(B-13) 


(B-14) 


Assuming the solutions of the form exp(-jp-r_) / the dis- 


persion equation becomes 

p^ - (a+c)p^ + (ac-bd) = 0 


(B-15) 


Then the wave equation describing the transverse components 


of em waves is given by 


(V^ + p?)E . = 0 . 
^ T ^ 2 .' Z1 


(B-16) 


The longitudinal field components inside the annular plasma 
(a<r<b) are 






im<j)-rz 


(B-17) 
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and 

H = 
Z 


h,E T + h„E „ 
1 zl 2 z2 




AhiJ^(Pir) + Ehj^Nj^(Pj^r) + Ch2J^(P2r) + DIi2Hj^(P2r)] 


jm(j)-rz 


(B-18) 


where A, B, C and D are arbitrary, complex constants. 


Using 




- 1 /r 

R/b 

-Kj^/FK^ 

S/b 

' 


r/b 

P/b 

0 

Q/b 


I 



'^T 





^2 '^T 

E T + 
zl 

P 

1 T 

^z2 




^z^"t 








(B-19) 


t 


the transverse E field component inside the annular plasma 
column (a<r<b) are 

^2Pl^m^Pl^^} + cf- ^lP2‘^m ^P2^^ 

+ d{- ^iN^(P2^^) + . (B-20) 
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hPA‘Pl=^)} + '^’{^'‘lJm(P2P> + hP2'';'P2P>} 

+ D{^«iNj^(P2r) + Lj^PjNj;^(P2r)j]e3’"'''-'’^ . 


(B-21) 


For the transverse H field components: 


= H' + y2Pl'^i;,<PlP>}+ ®{- ^^2’'m<PlP> 

+ y2Pl’'m<Pl'^>? + ^Vm<P2p> + yiP2^m<P2p>} 

+ d{- + yiPaH^fPsP)}]®^”"’’’'''’ , (b-22) 


“ H^y2-^m<PlP> + ^2Pl'^m'Pl’'>} ■" ®{^y2\<PlP> 
^2Pl®m<PlP>i + ={“yiJm<P2^> + <P2P>) 

■" °i%A'P2P> + ^lP2«m<P2'^>}]®^”'*'''^ . 


(B-23) 


Now the boundary conditions at the plasma surface (r=a) 
shall be taken as the continuity of the tangential E and H 
fields; 

= UJj^(qa) , (B-24) 

+ ^Vm(Pl^) + Dh2Nni(P2^^ =VJ^(qa), 


(B-25) 
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- + V?Z^j;(qa) , 


+ I*2PlJm'Pl®>] + ‘Pl®> + ^2 PA'Pi^>} 

+ cl^^Vm'P2^> + ^lP2'=';<P2^)|+ + ^lP2''i'P2^>} 

(B-26) 

'^{^y2''m‘Pl»> + ®lf^y2\<Pl^> ■" ^2Pl<<Pl^>} 

+ '=(i^yi^i»<P2“> + ^lP2'^m<P2®>} + Pff^yA<P2^> •" ^iP2VP2»'} 

= - U?Y^J^{qa) + viH^J^(qa) . (B-27) 

Also at the boundary of the perfect conductor wavequide 
wall (r=b) , 

AJi^(Plb) + BN^(p^b) + CJ^(p2b) + DN^(p2b) = 0 , (B-28) 

+ _+ Cff Vm‘P2>^> P'lP2'"i<P2‘>>} 

+ P{^^l*'m<P2*>> ^ V2«m<P2'=>} = ° ' 

In order to have non-trivial solutions of A, B, C, D, 

U and V, the 6x6 determinant of these coefficients must 


be zero; 
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SjL (Pib) 

S^(p^a) 


S2(Pj^b) 

N^(P].a) 

S2(Pia) 

T2(Pj^a) 


■^m<P2*=> 
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T 3 (p^a) 


*P2*=.> 

S4 (P2b) 

•’m‘P2®’ 

Vm'P2®> 
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T4 (P2a) 


0 

0 

-Ji„(qa) 

0 


-^?J^(qa) 

;Y J' (qa) 
o m ^ 


0 

0 

0 

-J^(qa) 
-%J (qa) 


0 


. (B-30) 



APPENDIX C. MATCHING FIELDS AT AIR -PLASMA INTERFACE 


Let us consider the configuration shown in Fig. (3-1) 
when the incident wave is TE^^. 

Region 1 

All components of E field in region 1 are 


M 

M n 


E = 0 + I-R , J (q r)e 
z nl m^^n 

n 


(C-1) 


E 


jm(}i-r^z r 

= C Z J* Cq, r)e + I" 

o m ^1 - h 

n>- 


"R T [•=^5J (q r) j 
nl r 


jitu|)+r^z jm(}i+r^z 

•e ^ + [R^-,?Z J'(q r)]e ^ 

nl^ o m ^n 


. (C-2) 


and 


E 


■ „o^,Ev-'^ 1 ,v 

J g J (q,r)e +■ ) 

E m "1 




jiti(})+rj^z 

•e ^ - [R:,jC-^^(q:r)]e 


E ,.E ,, j-Of+A'’ 

■nP'-r^m‘V ■ 


q ^ 


CC-3) 


Similarly all components of H field in region 1 is 


.E. 


H = 


jm(|)-r^z 

J (q r)e + T-R -,J (q r) e 

m ^1 ^ nl m ^n 

n 




, (C-4) 



144 


H, = (q?r)e 

4> ‘ 




n 


[-CY J' Cg^r) ] 
nl o m 


;]itt(f>+r^z . jin^+r^z 

> +[-R (q r)]e 

nl r m 


, (C-5) 


and 


■ A thE 

■p p D^t ^ t ^ 

n 


. Y m ... 

r • r O T , M . , 

R , [D 5 - 73 — j (q r) ] 
nl M m ^n 
q r , 


M In' *1 

jm(!)-r z „ jitn|)+r z^ 

e + [-R^,q^?J' (q®r)]e 

^ nl^n^ m 

.E 


. (C-6) 


Here the propagation constant, is determined from 

the characteristic equation of TE wave, i.e., 


J' (q^b) = 0 
m ^n 


(C-7) 


-where 


2 2 

E ^ ^ 2 

^n = ^n + 


(C-8) 


,M , 


Also the propagation constant^ is determined from 
the characteristic equation of TM waves, i.e., 

= 0 , (C^9) 


when 


+ k 2 

^n n o 


(C-10) 
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Region 2 

A: In the plasma a<r<b, all components of F field 

in the plasma region are 


E = y T „ 

z ^ n2 
n 


+ 


^ jm(|)~r z P 

D N (p„r) e ^ + V R _ A J (FtT) 

n m ^2 M ^ n2j_ n m ^1 


jm<j)+r z. 


(C-11) 


= I T ^Fa {iS- J (p,r) + pTJ'(pTr)} 
? “ n2| n r 2 m^’‘=^l 2 1 m'^1 ' 


+ Pj_ N_;^(Pj_r)} 


+ Dj, {3f f-l N^tPjr) + P2 n; (p^r)}- 


n t. 


+ ^2PA <Pi=^>> + <‘^'"m<P2^> + 

P2'^A<P2=^>> V^^X<P2-^> ^!p2'^;‘P2^>> 


jm(f)+r^z / 


(C-12) 
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and 


E 


j '"n2[V-^^2 Pi ^m'PlP' 

n L 


+ ^2 '’m CPl"^' + "2 Pi (Pl"^’ 


+ C^t-iS lJ J^CPir) + J^^PjJ^CpjD) 


^1 ’'m‘P2Pl + ‘‘l P2®';'P2P>^ 


jm(f)-r z r ' -irr, 

P " +• I ^2 V-^ ^„'PlP> 

n . L 


+ ^2 Pp '^m*PlP' ■*■ ®n^“^r ^2 '*'m*PlP* 

+ '‘2 Pi ^1 "'m‘P2P> 


+ J.P pj j;CP2r)) + D„{-% LiJ' H„(P,r) 


n 


+ Up P 2 «;<?2p> 


jmd)+r z 
e-^ ^ n 


(C-13) 


Similarly all components of H field in the plasma 
region are 
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-, jm<})-r z 

+ Wm<P2>^i VX'P2=^>J® 

-, jm'4>+r z 

+ C h_J (p^r) + D h^N (p^r) e ^ / (C-14) 

n 2 m ^2 n 2 m ^2 J 

^ ^2 + ^2 . 
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■ + Jm‘P2>^' P2 ‘P2>^>) 
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\ = I ’’n2[V- ^ ""a ^2 % ’ 

n ^ 

+ E^{- 3| y\ H^CPj^r) + Y2 Pi N^tPl>^)> 

+ Or^{- 3B yt j^Cp^r, + yP p^ JitPa-^)} 
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E, = 


. jm^-r z 

^ T „ [U Cq r) + V ?Z J' (q r) ]e ^ 

^ n2 n r n . o m 
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+. y R „[-U Cq r) +■ V J' (q r)]e 
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IP. -n 

Similarly all components of H field in the center section 
are 


H = 


H. = 


-n- 
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jm<j)-r z jm<j)+r z 

''naVm^Vie ® 

(C-20-) 
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r ^ n2 ^ n-^ q r m ^n n^n^ m ^n 
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Y m jm<i)+r z 

+ y R [U jC — J (q r) - V q EJ'(q r) ] e ^ . 

(C-22) 


^ n n-!’q^r m'^n ‘ ' n.-‘n'=’~m ' ^n’ 


Here the propagation constant, F^, is determined from 
the zeros of 6x6 determinant obtained from the boundary 
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conditions , which is given in Appendix B . 


Region 3 

All components of E fields in the region 3 are 


M 

M ^ n 


E = I (q^\)e 

z ^ n3 m ^n 
n 


(C-23) 


Ea = 


. M • j t^E 
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Similarly all components of H fields in the region 3 are 


„ „ jmct)~r^z 


H = y T -J (q r)e 
z ^ n3 m ^n 
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(q^r)l 
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(C-27) 


and 
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• 1 tiE 

3 iti(j)-r z- 


Tangential components of E and H fields are matched at the 
air-plasma interface z=0, The reflection and trans- 

mission coefficients R's and T’s are determined by the 
point -matching method. 
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